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Abstract 

Unusual properties of strongly correlated liquid observed in the high-T c superconductors and heavy- 
fermion (HF) metals are determined by quantum phase transitions taking place at their critical points. 
Therefore, direct experimental studies of these transitions and critical points are of crucial importance for 
understanding the physics of high-T c superconductors and HF metals. In case of high-T c superconductors 
such direct experimental studies are absent since at low temperatures corresponding critical points are oc- 
cupied by the superconductivity. Recent experimental data on the behavior of HF metals illuminate both 
the nature of these critical points and the nature of the phase transitions. We show that it is of crucial 
importance to simultaneously carry out studies of both the high-T c superconductivity and the anomalous 
behavior of HF metals. The understanding of this fact has been problematic largely because of the absence 
of theoretical guidance. The main features of the fermion condensation quantum phase transition (FCQPT), 
which are distinctive in several aspects from that of conventional quantum phase transition (CQPT), are 
considered. Our paper deals with these fundamental problems through studies of the behavior of quasipar- 
ticles, leading to good quantitative agreement with experimental facts. We show that in contrast to CQPT, 
whose physics in the critical region is dominated by thermal and quantum fluctuations and characterized by 
the absence of quasiparticles, the physics of a Fermi system near FCQPT or undergone FCQPT is controlled 
by the system of Landau-type quasiparticles. However, contrary to the conventional Landau quasiparticles, 
the effective mass of these strongly depends on the temperature T, magnetic fields B, the number density x, 
etc. Our general consideration suggests that FCQPT and the emergence of novel quasiparticles near and be- 
hind FCQPT are distinctive features of strongly correlated substances such as the high-T c superconductors 
and HF metals. We show that the main properties and universal behavior of the high-T c superconductors 
and HF metals can be understood within the framework of presented here theory based on FCQPT. A 
large number of the experimental evidences in favor of the existence of FCQPT in high-T c superconductors 
and HF metals is presented. We demonstrate that the essence of strongly correlated electron liquids can 
be controlled by both magnetic field B and temperature T. Thus, the main properties of heavy-fermion 
metal such as magnetoresistance, resistivity, specific heat, magnetization, volume thermal expansion, etc, 
are determined by its position on the B — T phase diagram. The obtained results are in good agreement 
with recent facts and observations. 

PACS numbers: 71.27.+a, 74.20.Fg, 74.25. Jb 
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1. INTRODUCTION 



In the last two decades a new class of materials such as heavy-fermion (HF) metals and high-T c 
superconductors has been found, which display a dazzling variety of physical phenomena, see e.g. 

HQ. 

The high-T c superconductors are strongly correlated metals with normal state properties 
that are not at all those of a normal Landau Fermi liquid (LFL). In the case of HF metals, the 
electronic strong correlations result in a renormalization of the effective mass of the quasiparticles, 
which can exceed the bare mass by a factor up to 1000 or even diverge. These non-Fermi-liquid (NFL) 
systems demonstrate anomalous pure power-law temperature dependences in their low temperature 
properties over broad temperature ranges. It is believed that in this class the basic assumption of 
the LFL theory that at low energies the electrons in a metal should behave as weakly interacting 
quasiparticles is violated. Therefore, it is generally accepted that the fundamental physics that gives 
rise to the high-T c superconductivity and NFL behavior with a recovery of the LFL behavior under 
the application of magnetic fields observed in HF metals and high-T c compounds is controlled by 
quantum phase transitions. This has made quantum phase transitions a subject of intense current 
interest, see e.g. P, [3,0, [J- 

A quantum phase transition is driven by control parameters such as the composition, number 
density x or magnetic fields B, and takes place at a quantum critical point (QCP) when temperature 
T = 0. QCP separates an ordered phase generated by quantum phase transition from a disordered 
phase. It is expected that the universal behavior is only observable if the system in question is very 
near QCP, for example, when the correlation length is much larger than microscopic length scales. 
Quantum phase transitions of this sort are quite common, and we shall call them as conventional 
quantum phase transitions (CQPT). In the case of CQPT, the physics is dominated by thermal and 
quantum fluctuations of the critical state, which is characterized by the absence of quasiparticles. It 
is believed that the absence of quasiparticle-like excitations is the main cause of the NFL behavior 
and other types of critical behavior in the quantum critical region. On the base of scaling related 
to the divergence of the correlation length, one can construct the critical contribution to the free 
energy and evaluate the corresponding properties such as critical exponents, the NFL behavior, etc. 
P U 0, Is]. However along this way one may expect difficulties. For example, having the only 
critical contribution, one has to describe different types of the behavior exhibited by different HF 
■netals, see e,. QElQ. Note that HF ,neta ls a r e three-di— , S ee e,. flflQ) 

and, thus, the type of behavior cannot be related to the dimension. The critical behavior observed in 
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measurements on HF metals takes place up to rather high temperatures comparable with the effective 
Fermi temperature Tj,. For example, the thermal expansion coefficient ot(T) measured on CeNi 2 Ge2 
shows a 1/VT divergence over more than two orders of magnitude in temperature drop from 6 K 
down to at least 50 mK jfj. It is hardly possible to understand such a behavior on the base of the 
assumption of scaling when the correlation length has to be much larger than microscopic length 
scales. Obviously, such a situation can take place only at T — > 0. At some temperature T ~ T^, this 
macroscopically large correlation length must be destroyed by thermal fluctuations. 

The next problem is related to explanations of the recovery of the LFL behavior under applied 



magnetic fields B, observed in HF metals and the high-T c compounds, see e.g. [l|, [10j, • At T — > 0, 
the magnetic field dependence of the coefficient A(B), causing an electrical resistivity contribution 
Ap = A(B)T 2 , the Sommerfeld coefficient 7(-B) and x{B) in specific heat, C/T — 7(-B), and magnetic 
susceptibility, x(B), shows that A(B) ~ 1 2 {B) and A(B) ~ y 2 (B), so that the Kadowaki- Woods 
ratio, K = A(B)/ / y 2 (B) [9], is ^-independent and conserved lOj. Such a universal behavior is hardly 



possible to explain within the picture assuming the absence of quasiparticles which takes place near 
QCP of the corresponding CQPT. As a consequence, for example, these facts are in variance to 
;he spin-density- wave scenario ^ and the renormalization group treatment of quantum criticality 



12[. Moreover, striking recent measurements of the specific heat, charge and heat transport and the 



resistivity used to study the nature of magnetic field-induced QCP in heavy-fermion metal CeCoIn 5 



13 



14l | certainly seem to disagree with descriptions based on CQPT. 



In a system of interacting bosons at temperatures lower than the temperature of Bose-Einstein 
condensation a finite number of particles is concentrated in the lowest level. In the case of a 

noninteracting Bose gas at zero temperature, T = 0, this number is simply equal to the total number 
of particles in the system. In a homogeneous system of noninteracting Bosons, the lowest level is the 
state with zero momentum, and the ground state energy is equal to zero. For a noninteracting Fermi 
system such a state is impossible, and its ground state energy E gs reduces to the kinetic energy 
and is proportional to the total number of particles. Imagine an interacting system of fermions 
with a pure repulsive interaction. Let us increase its interaction strength. As soon as it becomes 
sufficiently large and the potential energy starts to prevail over the kinetic energy, we can expect the 
system to undergo a phase transition at T = 0. This takes place when the density x tends to the 
critical quantum point x — > Xpc an d the kinetic energy becomes frustrated, while the effective 
inter-electron interaction, or the Landau amplitude, being sufficiently large, starts to determine the 
occupation numbers of quasiparticles n(p) which deliver the minimum value to the ground state 
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energy E[n(p)]. As a result, at x < Xpc, the function n(p) is given by the standard equation that 
determines the minimum of functional J5[n(p)] 13, LLSl, uM I20L \21\ 

on(pj 

Here we deal with three dimensional (3D) case and assume that the phase transition takes place at 
x < xfc- At T = Eq. (1.1) determines the quasiparticle distribution function n (p), which delivers 
the minimum value to the ground state energy E. The function nn(p) being the signature of the 
new state of quantum liquids [22j does not coincide with the quasiparticle distribution function of 
the LFL theory in the region (pf — pi), so that < nn(p) < 1 and pi < pf < Pf, with pp = (S^x) 1 / 3 
being the Fermi momentum. Such a state was called the state with fermion condensate (FC) because 
quasiparticles located in the region (pf—pi) of momentum space are pinned to the chemical potential 
\x |itL IiqI 1^. We note that the behavior obtained for the single- par ticle spectrum and quasiparticle 
distribution functions is observed within exactly solvable models [23. 24, 25]. Lowering the potential 
energy, FC decreases the total energy. Unlike the Bose-Einstein condensation, which occurs even in 
a system of noninteracting bosons, FC can take place if the coupling constant of the interaction is 
large, or the corresponding Landau amplitudes are large and repulsive. 

We note the remarkable peculiarity of FCQPT at T = 0: this transition is related to the 
spontaneous breaking of gauge symmetry, when the superconducting order parameter n(p) = 
• (1 — uq(p)) has a nonzero value over the region occupied by FC, with the entropy 5 = 
26|| . At small values of the pairing coupling constant An, the gap A(p) is linear in An and 



vanishes provided that An — > 0, while «(p) remains finite 20j, |2lJ. As we shall see, this peculiarity 
allows to construct the theory of high-T c superconductivity based on FCQPT. Thus, the state with 
FC cannot exist at any finite temperatures and is driven by the parameter x: at x > xpe the system 
is on the disordered side of FCQPT; at x — xpci Eq. (1.1) possesses the non-trivial solutions nn(p) 
with pi = pf — Pf', while at x < xfc, the system is on the ordered side [26]. 

One of the most challenging problems of modern condensed matter physics is the structure and 
properties of Fermi systems with large coupling constants. The first solution to this problem was 
offered by the Landau theory of Fermi liquids, later called "normal", by introducing the notion of 
quasiparticles and so called amplitudes, which characterize the effective interaction among them j^. 
The Landau theory can be viewed as the low energy effective theory in which high energy degrees of 
freedom are removed by introducing the effective amplitudes instead of strong inter-particle interac- 
tion. Usually, it is assumed that the stability of the ground state of Landau liquid is determined by 
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the Pomeranchuk stability conditions: the stability is violated when even one of the Landau effective 
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28| . Note that the new phase, at 



interaction parameters is negative and reaches a critical value 
which the stability conditions are restored, can in principle be again described within the framework 
of the same theory. However, it has been demonstrated rather recently [ij 7 ] that the Pomeranchuk 
stability conditions cover not all possible instabilities: one of them is missed. It corresponds to the 
situation when, at the temperature T = 0, the effective mass, the most important characteristic 
of Landau quasiparticles, can become infinitely large. Such a situation, leading to profound con- 
sequences, can take place when the corresponding Landau amplitude being repulsive reaches some 
critical value. This leads to a completely new class of strongly correlated Fermi liquids with FC 



221 ]. which is separated from that of a normal Fermi liquid by the fermion condensation 



quantum phase transition (FCQPT) 



3Q 



In the FCQPT case we are dealing with the strong coupling limit where an absolutely reliable 
answer cannot be given on the bases of pure theoretical first principle foundation. Therefore, the 
only way to verify that FC occurs is to consider both exactly solvable models and experimental 
facts, which can be interpreted as confirming the existence of such a state. The exactly solvable 
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On the 



models unambiguously demonstrate that Fermi liquids with FC do exist, see e.g. [23J, 
other hand, these facts are seen in features of those two-dimensional (2D) systems with interacting 
electrons or holes, which can be represented by doped quantum wells and high-T c superconductors. 
Considering the HF metals and the 2D systems of 3 He, we will show that FC exist also in these 
systems. 

The goal of our review is to describe the behavior of Fermi systems with FC and to show that 
the existing data on strongly correlated liquids represented by the electronic (or hole) systems of 
high-T c superconductors and HF metals can be well understood within the theory of Fermi liquids 
based on FCQPT. In Section 2, we review the general features of Fermi liquids with FC in their 
normal state. Section 3 is devoted to consideration of the superconductivity in the presence of FC. 
We show that the superconducting state is totally transformed by the presence of FC. For instance, 
the maximum value Ai of the superconducting gap can be as large as Ai ~ O.Ief, while for normal 
superconductors one has Ai ~ 10 -3 £p, where Ep is the Fermi level energy. In Section 4 we describe 
the quasiparticle's dispersion and its lineshape and show that they strongly deviate from the case of 
normal Landau liquids. In Section 5 we consider the field-induced LFL in the heavy electron liquid 
with FC. In Section 6 we apply our theory to explain the main properties of magnetic-field induced 
Landau Fermi liquid in the high-T c metals. Section 7 is devoted to the appearance of FCQPT in 
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different Fermi liquids. In Section 8 we analyze the main properties of HF metals whose electronic 
system is placed on the disordered side of FCQPT. HF metals with the electronic system located on 
the ordered side of FCQPT are considered in Section 9. In Section 10 we show that FC manifests 
itself in the dissymmetry of tunnelling conductivity which can be observed in measurements on the 
high-T c compounds and HF metals. Finally, in Section 11, we summarize our main results. 



2. FERMI LIQUIDS WITH FERMION CONDENSATE 

To study the universal behavior of the high-T c superconductors and HF metals at low tempera- 
tures, we use the heavy electron liquid model in order to ignore the complications of the anisotropy 
of the lattice of solids and its microscopic inhomogeneity. It is possible since we consider the univer- 
sal behavior demonstrated by these materials and processes related to the power-low divergences of 
observables such as the effective mass, specific heat, thermal expansion, etc. These divergences are 
determined by small momenta transferred as compared to momenta of the same order of magnitude 
as those of the reciprocal lattice cell, and contributions coming from them can be safely ignored. On 
the other side, we can simply use the common concept of the applicability of the LFL theory when 



describing electronic properties of metals |27j. Thus, we may usefully ignore the complications due 
to lattice and its anisotropy. As a result, we regard the medium as homogeneous heavy electron 
isotropic liquid. 



2.1. Landau Theory of Fermi Liquid 



Let us start by explaining the important points of the LFL theory |27j. The LFL theory rests on 
the notion of quasiparticles which represent elementary excitations of a Fermi liquid. Therefore these 
are appropriate excitations to describe the low temperature thermodynamic properties. In the case 
of an electron system, these are characterized by the electron's quantum numbers and effective mass 
M*. The ground state energy of the system in question is a functional of the quasiparticle occupation 
numbers (or quasiparticle distribution function) n(p, T), just like the free energy F[n(p, T)], entropy 
S[n(p, T)], and other thermodynamic functions. From the condition that the free energy F = E — TS 
should be minimal, we can find the distribution function 

$(F - fiN) , T , m T] l-n(p,T) 

— — — = e p,T )-n T -Tin — - — = 0. 2.1 

5n(p,T) n(p,T) 
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Here \x is the chemical potential, while 

, T) = SE[n( P ,T)] 

lP ' ' 6n(p,T) ' 1 ] 

is the quasiparticle energy. This energy is a functional of n( p, T ) just like the total energy E[n(p, T)]. 
The entropy S[n(p,T)] is given by the familiar expression |2j| 

S[n{p, T)] = -2 J [n(p, T) In n{p, T) + (1 - n(p, T)) ln(l - n(p, T))} -0-, (2.3) 

which stems from purely combinatorial considerations. Equation (2.1) is usually presented as the 



Fermi-Dirac distribution 

-i 



n(p,T) = <l + exp 



(g(p,T)-/x) 
T 



(2.4) 



At T — > 0, one gets from Eqs. (2.1) and (2.4) the standard solution n(p,T —>■ 0) — > #(pf — p), with 
6*(pir — p) is the step function, e(p ~ p F ) — \i = Pf(p — Pf)I^li where M£ is the Landau effective 
mass |27| 

J_ _ 1 &(p,T = 0) , 

M£ ~ p dp lp=PF ' { ' 

It is implied that M£ is positive and finite at the Fermi momentum pp. As a result, the T-dependent 
corrections to M£, to the quasiparticle energy e(p), and to other quantities, start with T 2 -terms. 
The effective mass is given by the well-known Landau equation 

Here F a>ai (p-p, Pi) is the Landau amplitude depending on the momenta p, spins cr, and M is the bare 
mass of an electron. For the sake of simplicity, we omit the spin dependence of the effective mass 
since in the case of a homogeneous liquid and weak magnetic fields M£ does not noticeably depend 
on the spins. 

Applying Eq. (2.6) at T = and taking into account that n(p,T = 0) becomes the step function 
9(pf — p), we obtain the standard result 

Ml 1 



Here N is the density of states of the free Fermi gas and F 1 (pp 1 pp) is the p-wave component of 
the Landau interaction amplitude. Since in the LFL theory x = p F /37r 2 , the Landau amplitude 
can be written as F 1 (pf,Pf) = F l (x). Assume that at some critical point xfc the denominator 
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1 — N F 1 (pf,Pf)/3) tends to zero, that is (1 — N Q F 1 (x)/3) oc (x — xfc) + a (x — %fc) 2 + 



As a result, one obtains that M£(x) behaves as 



31 



32] 



M 



A 



x 



B 1 
oc -. 

- %fc r 



(2.7) 



Here A and B are constants and r = {x — xfc*) is the "distance" from the QCP of FCQPT taking 
place at xfc ,- The observed behavior is in g ood agreement with recent experimental observations (see 



e.g. 



33l l34j). and calculations 
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36l |37j, see Section 7 as well. In the case of electronic systems, 



Eq. (2.7) is valid at x > xfc when r > [19J, |38j. Such a behavior of the effective mass can be 
observed in the HF metals with a quite flat, narrow conduction band, corresponding to the large 
effective mass M* L [x ~ xfc\ with strong electron correlations and the effective Fermi temperature 



Tfc ~ p 2 F /Ml(x) of the order of a few Kelvin or even lower (see e.g. 



)■ 



2.2. Fermion condensation quantum phase transition 

As we have seen above at T = when r = (x — xpc) ~~ > 0, the effective mass diverges, M£(r) — > oc, 
and eventually beyond the critical point xpc the distance r becomes negative making the effective 
mass negative, as it follows from Eq. (2.7). To escape the possibility of being in unstable and in 
essence meaningless states with negative values of the effective mass, the system is to undergo a 
quantum phase transition at the critical point x = xpc- Because the kinetic energy near the Fermi 



surface is proportional to the inverse effective mass, t 
kinetic energy and can be recognized as FCQPT (29 



lis phase transition is triggered by the frustrated 



391 ] . Therefore behind the critical point xpc of 



this transition, the quasiparticle distribution represented by the step function does not minimize the 
Landau functional J5[n(p)]. As a result, at x < xfc the quasiparticle distribution is determined by 
the standard equation that is used to search for the minimum of the energy functional 17 1 

SE\n(p)\ , \ , N 

^^L = ,(P) = «P,<P<P/. (2-8) 

Equation (2.8) determines the quasiparticle distribution function no(p), which minimizes the ground 
state energy E. Being determined by Eq. (2.8), no(p) does not coincide with the step function in the 
region (pf—pi), so that < n (p) < 1, while outside the region it coincides with the step function. It 
follows from Eq. (2.8) that the single particle spectrum or the band is completely flat over the region. 
Such a state was called the state with FC because quasiparticles located in the region (pf — pi) of 
momentum space are pinned to the chemical potential fi |l7l . Ilfll |22| . We note that this behavior was 
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obtained for the band and quasiparticle distribution functions using some exactly solvable models 

3, Ei. 



*(P) 






n o<P) 



Pi P F 



FIG. 1: The quasiparticle distribution function no(p) and energy e(p). Since no(p) is the solution of Eq. 
(2.8) it implies no(p < pi) = 1, no(j»i < p < Pf) < 1 and no(p > p/) = 0, while e(pi < p < pj) = fx. The 
Fermi momentum pp obeys the condition p\ < pp < pf. 

The possible solution no(p) of Eq. (2.8) and the corresponding single particle spectrum e(p) are 
shown in Fig. 1. 

As we shall see in Section 3 the relevant order parameter of the FC state, *c(p) = 



yJio(p)(l ~ n o(p)), is the order parameter of the superconducting state with the infinitely small 
value of the superconducting gap. Therefore the entropy of this state is zero, S(T = 0) = 
Thus, this state is of pure quantum nature and cannot exist at any finite temperatures. This quan- 
tum state with FC is driven by the density x: at x > xfc the system is on the disordered side of 
FCQPT; at x = Xpc, Eq. (2.8) has the non-trivial solutions n (p) with p { = p F = pj; at x < Xpci 
the system is on the ordered side. We note that the solutions no(p) of Eq. (2.8) can be viewed as 
new to the LFL theory. Indeed, at T = 0, the standard solution n(p, T —>■ 0) —>■ 6(pf — V) is n °f 
the only one possible. The "anomalous" solutions of Eq. (2.1) can exist because the logarithm on 
the right hand side of Eq. (2.1) is finite when p belongs to the region (pj — p^ and < n (p) < 1. 
Therefore in the region Tln[(l — uq(p, T))/no(p, T)]\t->o ~^ 0, an d we again arrive at Eq. (2.8). 

Let us assume that with the decrease of the density (or with the growth of the interaction strength) 
FC has just taken place. It means that Pi — > pj — > pp, and the deviation Sn(p) = n (p) — 9(pp — p) 
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is small. Expanding the functional i£[n(p)] in Taylor's series with respect to Sn(jp) and retaining the 
leading terms, one obtains from Eq. (2.8) the following relation 



e(p) = £ (p) + / F(p, Pl )5n( Pl )^ ; p t < p < Pf , (2.9) 



where F(p, px) = 5 2 E/5n(p)5n(pi) is the Landau amplitude. Both quantities, the interaction and 
the single-particle energy £o(p) are calculated at n(p) = 0{j>F—p). Equation (2.9) acquires nontrivial 
solutions at some density x = xpc- Thus, FCQPT takes place if the Landau amplitudes depending 
on the density are positive and sufficiently large, so that the potential energy is bigger than the 
kinetic energy. Then the transformation of the Fermi s tep function n(p) = 6(pf — p) into the 



smooth function defined by Eq. (2.9) becomes possible |17|, The system in question can be 
considered as a strongly correlated Fermi liquid behind the point Xpc- It is seen from Eq. (2.9) that 
the FC quasiparticles form a collective state, since their energies are defined by the macroscopical 
number of quasiparticles within the momentum region (pf—pi)- The shape of the excitation spectra 
related to FC is not affected by the Landau interaction, which, generally speaking, depends on the 
system's properties, including the collective states, the irregularity of the composition, impurities, 
etc. The only thing determined by the interaction is the width of the FC region (p/ — P i) provided 
the interaction is sufficiently strong to produce the FC phase transition at all. Thus, we can conclude 
that the spectra related to FC are of a universal form, being dependent, as we shall see in Subsections 
2.3 and 3.1, on the temperature and the superconducting gap. The existence of such spectra can be 



viewed as the characteristic feature of the "quantum protectorate" 



40 



4l|. 



2.3. The "shadow" of the fermion condensate at finite temperatures 

According to Eq. (2.1), the single-particle energy e(p, T) within the interval (pf — Pi) at T <C Tf 



is linear in T 



421 ] . At the Fermi level, one obtains by expanding ln(...) in terms of n(p) 



e(p,T)-^(T) = Tin — — ~ T — — . 2.10 

n(p) n(p) 



Here Tf is the temperature, above which FC effects become insignificant 

Z> „ ±^L „ 5££ . (2 . n) 

e F 2Me F VL F v 1 

In this formula VLpc is the FC volume, Bp is the Fermi energy, and VLp is the volume of the Fermi 
sphere. We note that at T Tf the occupation numbers n(p) are approximately independent of T, 
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being given by Eq. (2.8). At finite temperatures according to Eq. (2.10), the dispersionless plateau 
e(p) = \i is slightly turned counter-clockwise about /i. As a result, the plateau is just a little tilted 
and rounded off at the end points. According to Eqs. (2.5) and (2.10), the effective mass M FC 
related to FC is given by, 

M* FC * Pf PJ ^ ■ (2.12) 

To obtain Eq. (2.12) an approximation for the derivative dn(p)/dp ~ —l/(pf — Pi) was used. It is 
seen from Eq. (2.12) that at < T Tf, the heavy electron liquid behaves as if it were placed 
at QCP, in fact it is placed at the quantum critical line x < xfc-, that is the critical behavior is 
observed at T — > for all x < xfc- 

Having in mind that (pf — pi) <C Pf and using Eqs. (2.11) and (2.12), the following estimates for 



the effective mass M FC are obtained: 



M* FC N(0) T f 



(2.13) 



M iV (0) T 

Eqs. (2.12) and (2.13) show the temperature dependence of M FC . In Eq. (2.13) iV (0) denotes the 
density of states of noninteracting electron gas, and N(0) is the density of states at the Fermi level. 
Multiplying both sides of Eq. (2.12) by (pf —Pi), we obtain the energy scale E separating the slow 
dispersing low energy part related to the effective mass MZp from the faster dispersing relatively 
high energy part defined by the effective mass M£ (29, 30 



E ~ 4T . (2.14) 

It is seen from Eq. (2.14) that the scale E does not depend on the condensate volume. The single 
particle excitations are defined according to Eq. (2.12) by the temperature and by (pf — pi), given 
by Eq. (2.8). Thus, we conclude that the one-electron spectrum is of universal form and has the 
features of the " quantum protectorate" . 

It is pertinent to note that outside the FC region the single particle spectrum is not strongly 
affected by the temperature, being defined by M£. Thus, we come to the conclusion that a system 
with FC is characterized by two effective masses: M FC which is related to the single particle spectrum 
at lower energy scale and M£ describing the spectrum at higher energy scale. The existence of 
two effective masses is manifested by a break (or kink) in the quasiparticle dispersion, which can 
be approximated by two straight lines intersecting at the energy Eq. This break takes place at 



temperatures T c < T <C Tf, which is in accord with experimental data Q], and, as we will see, 
at T < T c which is also in accord with the experimental facts 0, |^. Here T c is the critical 
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tempereture of the superconducting phase transition. The quasiparticle formalism is applicable to 
this problem since the width 7 of single particle excitations is not large compared to their energy, 
being proportional to the temperature, 7 ~ T at T > T c j^. The lineshape can be approximated 



by a simple Lorentzian 43, . This is consistent with experimental data obtained from scans at a 
constant binding energy |46( (see Sec. 4). 

It is essential to have in mind, that the onset of the charge density wave instability in a many- 
electron system, such as an electron liquid, which takes place as soon as the effective inter-electron 
constant reaches its critical value r s = r cdw , is preceded by the unlimited growth of the effective 
mass, which is demonstrated Section 7. Here r s = r^jas with r$ being the average distance between 
electrons, while a# is the Bohr radius. Therefore the FC occurs before the onset of the charge density 
wave. Hence, at T = 0, when r s reaches its critical value r FC corresponding to Xpci r FC < r cdw , 



FCQPT already inevitably takes place 



381 ] . It is pertinent to note that this growth of the effective 



mass with decreasing electron density was observed experimentally in a metallic 2D electron system 
in silicon at r s ~ 7.5 [s^j]. Therefore we can take this value as an estimate, r F o ~ 7.5. On the other 
hand, there exist ch arge density waves or strong fluctuations of charge ordering in underdoped high- 
T c superconductors j^j , . Thus, the formation of FC in high-T c compounds can be thought more 
as a general property of low density electron liquid embedded in these solids rather than an unusual 
and anomalous solution of Eq. (2.8) |38j]. Beyond the point of FCQPT, the condensate volume is 
proportional to {r s — r FC ) as well as Tf/sp ~ (r s — r FC )/r FC at least when {r s — r FC )/r FC <C 1. 
Therefore we obtain 

r s -r FC ^ p f - Pl ^ x FC -x _ 
r FC Pf Xfc 

FC serves as a stimulator that creates new phase transitions, which eliminates the degeneracy of the 

spectrum. For example FC can generate spin density waves or antiferromagnetic phase transition, 

thus leading to a whole variety of new properties of the system under consideration. Then, the 

onset of the charge density wave is preceded by FCQPT, and both of these phases can coexist at 

the sufficiently low density when r s > r c d w . The transition to superconductivity is strongly assisted 

by FC because both of the phases are characterized by the same order parameter. As a result, the 

superconductivity by removing the spectrum degeneracy "wins" the competition with other phase 

transitions up to the critical temperature T c , see Section 3. We now turn to the consideration of the 

superconducting state and quasiparticle dispersions at T < T c . 
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3. THE SUPERCONDUCTING STATE WITH FERMION CONDENSATE 

In this Section we consider the superconducting state of 2D heavy electron liquid since the high-T c 
superconductors are predominantly represented by 2D structures. On the other hand, our consider- 
ation can be easily adopted to the 3D case. 



3.1. Superconducting state at T = 

At T = 0, the ground state energy E 9S [k(p), n(p)] of a 2D electron liquid is a functional of both 
the order parameter of the superconducting state n(p) and the quasiparticle occupation numbers 
n(p). This energy is determined by the well-known equation of the Bardeen-Cooper-Schrieffer (BCS) 



weak-coupling theory of superconductivity, see e.g. |4£ 



. 50] 



E gs = E[n(p)} + A / n Pl ,p 2 )K( Pl K(p 2 )^^ • (3.1) 
Here ^^(p)] is the ground-state energy of a normal Fermi liquid, and 



n(p) = v (p); re(p) = v(p)u(p) = yn(p)(l - n(p)), (3.2) 

Where u(p) and v(p) are the normalized coherence factors, v 2 (p) + u 2 (p) = 1. It is assumed that 
the pairing interaction AoV(pi,P2) is weak. We define the superconducting gap 

A(p) = -A / V(p,ih)K(pi)^. (3.3) 

Minimizing E gs with respect to v(p) we obtain the equation connecting the single-particle energy 
e(p) to A(p), 

1 - 2v 2 (r>) 

e(p)-v = A(p) Pj , (3.4) 

Z/ HPJ 

where the single-particle energy e(p) is determined by the Landau equation (2.2). With insertion of 
the value of n(p) into Eq. (3.3), we obtain the known equation for A(p) 

A (p)= - T ynp.P0 7 ^ n?7 ^ ^. (3.5) 

If Ao — > 0, then the maximum value Ai of the superconducting gap A(p) tends to zero, and Eq. 
(3.4) reduces to Eq. (2.8) 

5Mn(p)l 

= e(p) — fj, — 0, provided that < n(p) < 1; or «(p) f 1 0; Pi < p < p f . (3.6) 

on^pj 
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It is seen from Eq. (3.6) that at x < xfc-, the function n(p) is determined by the standard equation to 
search the minimum „ t — E lnip)] 0Q. Elation ,3,, detunes the ^particle dis- 
tribution function no(p) which delivers the minimum value to the ground state energy £ , gs [ft(p), n(p)] 
in the Ao — > limit. Now we can study the relationships between the state defined by Eq. (3.6), or 
by Eq. (2.8), and the superconducting state. At T = 0, Eq. (3.6) defines the particular state of a 
Fermi-liquid with FC, for which the modulus of the order parameter |/«(p)| has finite values in the 
range of momenta Pi < p < Pf, while A x — > 0. Such a state can be considered as superconducting, 
with an infinitely small value of A 1; so that the entropy of this state is equal to zero. It is obvious 
that this state, being driven by the quantum phase transition, disappears at T > 0,13- Any 
quantum phase transition, which takes place at temperature T = 0, is determined by a control pa- 
rameter other than the temperature, for instance, by pressure, magnetic field, or the number density 
of mobile charge carriers x. In the case of FCQPT, as we have shown in Section 2, the control 
parameter is the density x of the system, which determines the strength of the Landau amplitudes. 
FCQPT occurs at the quantum critical point x = xfc- 

As any phase transition, FCQPT is related to the order parameter, which induces a broken 
symmetry. As it follows from Eq. (3.2), the order parameter of the state with FC is ft(p). Thus, 
the solutions n (p) of Eq. (3.6) represent a new class of the solutions of both BCS equations and 
LFL equations. In contrast to the conventional BCS solutions |49|, the new ones are characterized 
by the infinitesimal value of the superconducting gap, Ai — > 0, while the order parameter ft(p) is 
finite. At the same time, in contrast to the standard solutions of the LFL theory, the new ones 
are characterized by the superconducting order parameter n(p), that is at T — > 0, the quasiparticle 
distribution function does not tend to the step function 9(pf — p), being the solution of Eq. (3.6). 
Thus, we can conclude that the solutions of Eq. (3.6) can be viewed as common solutions of both 
BCS and LFL equations, while Eq. (3.6) can be derived starting from either BCS or LFL theory. As 
we shall see, on the basis of the peculiarities of this new class of solutions, we can define the notion of 
the strongly correlated Fermi liquid and explain the main properties of the high-T c superconductivity 
and HF metals. It is essential that the existence of state with FC can be revealed experimentally 
since the order parameter n(p) is suppressed by a magnetic field B. Destroying the state with FC, 
the magnetic field B converts the strongly correlated Fermi liquid into the normal LFL. In that case, 
the magnetic field plays the role of control parameter. 

When pi —>■ pf —>■ Pf, Eq. (3.6) determines the critical density xfc at which FCQPT takes place. 
When x < xfc, the system becomes divided into two quasiparticle subsystems: the first subsystem is 
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in the (pf —pi) range and is characterized by the quasiparticles with the effective mass M FC oc 1/Ai, 
as it follows from Eq. (3.4), while the second one is occupied by quasiparticles with finite mass M£ 
and momenta p < pi. When Ao — > 0, the density of states near the Fermi level tends to infinity, 
N(0) ~ M FC ~ 1/Ai. The quasiparticles with M FC occupy the same energy level and form pairs 
with binding energy of the order of Ai and with average momentum p , Po/pf ~ {Pf — Pi)/pf ^ 1- 
Therefore, this state strongly resembles the Bose-Einstein condensate, in which quasiparticles occupy 
the same energy level. But the FC quasiparticles have to be spread over the range (pf — pi) in 
momentum space due to the exclusion principle. In contrast to the Bose-Einstein condensation, the 
fermion condensation temperature is T c = 0. And in contrast to the ordinary superconductivity, 
FC is formed due to the Landau repulsive interaction F(p, pi) rather than by the relatively weak 
attractive quasiparticle-quasiparticle pairing interaction A ^(pi,P2)- 

If Ao 7^ 0, Ai becomes finite, leading to a finite value of the effective mass M FC , which can 
be obtained from Eq. (3.4) upon differentiating both parts of this equation with respect to the 
momentum p and using Eq. (2.5) [29! Isfli 

M FC * p F V -^ . (3.7) 

As to the energy scale, it is determined by the parameter Eq: 

E = e(p f ) - e( Pi ) ~ Pf&=^- * 2A 1 . (3.8) 

IVlpQ 

3.2. Superconducting state at finite temperatures 

Let us assume that the range of FC is small, that is (pj_— Pi)/PF ^ 1> an d 2Ai Tf so that 



the order parameter fc(p) is governed mainly by FC [29|, |30|. To solve Eq. (3.5) analytically, we 
take the BCS approximation for the interaction X V(p, pi) = — A if |^(p) — A*| < wc, i.e. the 
interaction is zero outside this region, with lud being the characteristic energy, e.g. that of phonon. 
As a result, the gap becomes dependent only on the temperature, A(p) = Ai(T), being independent 
of the momentum, and Eq. (3.5) takes the form 

E /2 



1 



N FC X I , ^ +N L X [ di (3.9) 



Here we set ^ = e(p) — \i and introduce the density of states N FC in the (p/ — pi) range, or E Q range. 
It follows from Eq. (3.7), that N F c = (j>f — p F )p F /2irAi(0). The density of states N L in the range 
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(uo D — Eq/2) has the standard form N L = Ml/2n. If the energy scale E — > 0, Eq. (3.9) reduces to 
the BCS equation. On the other hand, assuming that E < 2oj d and omitting the second integral 
on the right hand side of Eq. (3.9), we obtain 

Ax(0) = jgg%Zgj Infl + v^) = 2(3e F P -l^ In f 1 + V2) , (3.10) 

where the Fermi energy e F — pf,/2M£, and the dimensionless coupling constant (3 is given by the 
relation (3 = AoM£/27r. Taking the usual values of f3 for conventional superconductors as (3 ~ 0.3, 
and assuming (pf —Pf)/pf — 0.2, we get from Eq. (20) a large value of A^O) ~ O.lep, while for 
conventional superconductors one has a much smaller gap, Ax(0) ~ 10~ 3 ef- Taking into account the 
omitted above integral, we obtain 

M0) * 2^^^1n(l + V2)fl + /31n^). (3.11) 
Pf V Eq J 

It is seen from Eq. (3.11) that the correction due to the second integral is small, provided E ~ 2u D . 

Below we shall show that 2T C ~ A^O), which leads to the conclusion that the isotope effect is absent 

since Ai is independent of ujd- But this effect is restored as Eq — > 0. 

Assuming E ~ ud but E > cu D , we see that Eq. (3.11) has no standard solutions A(p) = 

Ai(T = 0) because ud < s(p — Pf) — H and the interaction vanishes at these momenta. The only 

way to obtain solutions is to restore the condition E < ujd- For instance, we can define such a 

momentum Pd < Pf that 

Ax(0) = 2(3e F PD ~ PF \n(l + V2) = u D , (3.12) 
p F v ' 

while the other part in the (j)f — Pi) range can be occupied by a gap A 2 of the different sign, 
Ai/A 2 < 0. It follows from Eq. (3.12) that the isotope effect is preserved, while both gaps can have 
s-wave symmetry. 

At T ~ T c , Eqs. (3.7) and (3.8) are replaced by Eqs. (2.12) and (2.14), which is valid also at 
T c < T « T f 

M FC ~ Vf VJ -^, E ~ 4T C ; if T c < Tthen , M FC ~ Pf^^, E ~ 4T . (3.13) 
Equation (3.9) is replaced by its conventional finite temperature generalization 

rE /2 % Je + A\{T) 



1 = N FC \o / , = = tanh V " — 
V / e + A?(T) 2T 



+ 



+ N L X / - I tanh V ' U . (3.14) 



/ Eo /2 ^2 + A 2 (T) 
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Putting Ai(T — > T c ) — > 0, we obtain from Eq. (3.14) 

2T C ~ Ai(0) , (3.15) 

with Ai(T = 0) being given by Eq. (3.11). Upon comparing Eqs. (3.7), (3.8) and (3.13), (3.15), we 
see that M FC and Eq are almost temperature independent at T < T c . 



3.3. Bogoliubov quasiparticles 

It is seen from Eq. (3.5) that the superconducting gap depends on the single-particle spectrum 
e(p). On the other hand, it follows from Eq. (3.4) that e(p) depends on A(p), provided that at 
Ao — > Eq. (3.6) has the solution determining the existence of FC. Let us assume that Ao is small 
enough, so that the particle-particle interaction XoV(p, pi) can only lead to a small perturbation of 
the order parameter n(p) determined by Eq. (3.6). It follows from Eq. (3.7) that the effective mass 
and the density of states N(0) oc M FO oc 1/Ai are finite. As a result, we conclude that in contrast 
to the conventional theory of superconductivity the single-particle spectrum e(p) strongly depends 
on the superconducting gap and we have to solve Eqs. (2.2) and (3.5) in a self-consistent way. On 
the other hand, let us assume that Eqs. (2.2) and (3.5) are solved, and the effective mass M FC is 
determined. Now one can fix the quasiparticle dispersion e(p) by choosing the effective mass M* 
of the system in question equal to the obtained M FC and then solve Eq. (3.5) as it is done in the 
case of the conventional BCS theory of superconductivity As a result, one observes that the 
superconducting state is characterized by the Bogoliubov quasiparticles (BQ) [51] with the dispersion 



£(p) = V( £ (p)-^) 2 + A2 (p)> 

and the normalization condition for the coherence factors v(p) and u(p) is held. We arrive to the 
conclusion that the observed features agree with the behavior of BQ predicted from the BCS theory. 
This observation suggests that the superconducting state with FC is BCS-like and implies the basic 
validity of the BCS formalism in describing the superconducting state 5^] . Although the maximum 
value of the superconducting gap given by Eq. (3.11) and other exotic properties are determined by 
the presence of the fermion condensate. It is exactly the case that was observed experimentally in 
the high-T c cuprate Bi 2 Sr2Ca2Cu30io+<5 j&J. 

We have returned back to the LFL theory since the high energy degrees of freedom are eliminated 
and the quasiparticles are introduced. The only difference between the Landau Fermi-liquid, which 
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serves as a basis when constructing the superconducting state, and Fermi liquid after FCQPT is that 
we have to expand the number of relevant low energy degrees of freedom by introducing a new type 
of quasiparticles with the effective mass M FC given by Eq. (3.7) and the energy scale Eq given by 
Eq. (3.8). Therefore, the dispersion e(p) is characterized by two effective masses M£ and M FC and 
by the scale En, which define the low temperature properties including the line shape of quasiparticle 

rsion of BQ has the standard form. We note that both the 



excitations \2{ 



effective mass M FC and the scale E Q are temperature independent at T < T c , where T c is the critical 
temperature of the superconducting phase transition 54J. At T > T c , the effective mass M FC and 
the scale Eq are given by Eqs. (2.12) and (2.14) respectively. Obviously, we cannot directly relate 
these new Landau Fermi-liquid quasiparticle excitations with the quasiparticle excitations of an ideal 
Fermi gas because the system in question has undergone FCQPT. Nonetheless, the main basis of the 
Landau Fermi liquid theory survives FCQPT: the notion of order parameter is preserved, and the 
low energy excitations of the strongly correlated liquid with FC are quasiparticles. 

As it was shown above, properties of these new quasiparticles are closely related to the properties 
of the superconducting state. We may say that the quasiparticle system in the range (pf—Pi) becomes 
very "soft" and is to be considered as a strongly correlated liquid. On the other hand, the system's 
properties and dynamics are dominated by a strong collective effect having its origin in its proximity 
to FCQPT and determined by the macroscopic number of quasiparticles forming FC in the range 
(j)f — Pi). Such a system cannot be perturbed by the scattering of individual quasiparticles and has 



features of a "quantum protectorate" and demonstrates the universal behavior |2j 
few remarks related to the quantum protectorate and the universal behavior j^J are in order here. 
As the Landau theory of Fermi liquid, the theory of the high-temperature superconductivity based 
on FCQPT deals with the quasiparticles which are elementary excitations of low energy. As a result, 
this theory gives general qualitative description of both the superconducting and normal states. Of 
course, one can choose the phenomenological parameters and obtain the quantitative consideration 
of the superconductivity as it can be done in the framework of the Landau theory when describing a 
normal Fermi-liquid, say 3 He. Therefore, we conclude that any theory, which is capable to describe 
FC and compatible with the BCS theory, will give the qualitative picture of the superconducting 
state and the normal state that coincides with the picture based on FCQPT. Both of the approaches 
can be agreed at a numerical level provided the corresponding parameters are adjusted. For example, 
since the formation of the flat band is possible in the Hubbard model [2j|, generally speaking one 
can repeat the results of the theory based on FCQPT in the Hubbard model. It is appropriate to 
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mention here that the correspondin g nu merical description limited to the case of T = has been 



obtained within the Hubbard model 



55 



56]. 



3.4. Pseudogap 

Now let us discuss some special features of the superconducting state with FC 0,0]. Consider 
two possible types of the superconducting gap A(p), namely that given by Eq. (3.5) and defined 
by the interaction XqV(p, pi). If this interaction is dominated by an attractive phonon-mediated 
attraction, the solution of Eq. (3.5) with the s-wave, or the s + d mixed waves will have the 
lowest energy. Provided the pairing interaction AoV(pi,p2) is the combination of both the attractive 
interaction and sufficiently strong repulsive interaction, the <i-wave superconductivity can take place 
see e.g. But both the s-wave symmetry and the ci-wave one lead to approximately the 

same value of the gap Ai in Eq. (3.11) 54]. Therefore the non-universal pairing symmetries in 
high-T c superconductivity are likely the result of the pairing interaction, while the (i-wave pairing 
symmetry is not essential. This point of view is supported by the data 0, 62, 63, 64, 65| . 

We can define the critical temperature T*, at which the superconductivity vanishes, as the tem- 
perature when A X (T*) =0. At T > T*, Eq. (3.14) has only the trivial solution A x = 0. On the other 
hand, the critical temperature T c can be defined as a temperature, at which the superconductivity 
disappears while the gap occupies only part of the Fermi surface. Thus, as we shall see there are 
;wo different temperatures T c and T* in the case of the d-wave symmetry of the gap. It was shown 



43, 



571 ] that in the presence of FC there exist nontrivial solutions of Eq. (3.14) at T c < T < T* when 
the BCS like interaction is replaced by the pairing interaction A ^(pi, P2), which includes the strong 
repulsive interaction leading to (i-wave superconductivity. In that case, the gap A(p) as a function 
of angle 4>, A(p) = AQop, (f>) possesses new nodes at T > T no d e , as it is illustrated by Fig. 2 0]. 

We show in Fig. 2 the ratio AQ»p, 4>)/T* calculated at three different temperatures 0.9T no d e , T no d e , 
and 1.2T no d e . An important difference between curve (a) and (b) and (c) is additional variations of 
the curves (b) and (c) marked by the two arrows. As seen from Fig. 2, the flattening occurs due to 
the appearance of the two new nodes. They appear at T = T n0( j, e and move apart with increase of the 
temperature confining the area 6 C shown by the two arrows in Fig. 2. It is also seen from Fig. 2 that 
the gap A is extremely small over the range 8 C . It was recently shown in a number of papers (see, 
e.g., 0,0]) that there exists an interplay between the magnetism and the superconductivity order 
parameters, leading to the damping of the magnetism order parameter below T c . And vise versa, 
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FIG. 2: The gap A(pp, <fi) as a function of (f> calculated at three different temperatures expressed in terms 
of T no de — T c , while A is presented in terms of T*. Curve (a), solid line, shows the gap calculated at 
temperature 0.9T„ Of fe. Curve (b), dashed line, presents the same at T = T no d e , and curve (c), dotted line, 
shows the gap calculated at temperature 1.2T no d e - Note the important difference in the curves (b) and (c) 
compared with curve (a) due to the flattening of the curves (b) and (c) about the nodes. The two arrows 
indicate the area 6 C emerging at T noc i e . 

one can anticipate the damping of the superconductivity order parameter by the magnetism. Thus, 
we conclude that the gap in the range 9 C can be destroyed by strong antiferromagnetic correlations 
(or by spin density waves) existing in optimally doped and underdoped superconductors. It is 
believed that impurities can easily destroy the gap A(p) in the considered area. As a result, the 
superconducting gap vanishes in the macroscopic area 9 C and causes the superconductivity to die out. 
We have to conclude that T c ~ T no d e , with the exact value of T c denned by the competition between 
the antiferromagnetic correlations (or spin density waves) and the superconducting correlations over 
the range 9 C . The behavior and the shape of the pseudogap is very similar to the ones of the 
superconducting gap as it is seen from Fig. 2. The main difference seen from Fig. 2 is that the 
pseudogap vanishes along segments 9 C of the Fermi surface, while the gap vanishes at the isolated 
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d-wave nodes. Our calculations show that the function 8 c (x) increases very fast at small values of x, 
9 c (x) ~ \/x. Therefore T c has to be close to T no d e - This result is in accord with numerical calculations 
of the function 8 C ([T — T c ]/T c ) plotted in Fig. 3. 



. 0,3 



FIG. 3: Calculated angle 6 C , pulling apart the two nodes, as a function of (T — T c )/T c . 

Thus, we observe that the pseudogap state appears at temperatures T > T c ~ T node and vanishes 
at T > T* when Eq. (3.14) has only the trivial solution Ai = 0. Quite naturally, one has to recognize 
that Ax scales with T*, while Eq. (3.15) takes the form 

2T* ~ Ai(0) . (3.16) 

It can then be concluded that the temperature T* has the physical meaning of the BCS transition 
temperature between the state with the order parameter k^O and the normal state. 

At T < T c quasiparticle excitations are characterized by sharp peaks. When temperature becomes 
T > T c and A(9) = in the range 9 C , there appear normal quasiparticle excitations with the width 7 
along the segments 9 C of the Fermi surface. There exists the pseudogap outside the segments 9 C , and 
the Fermi level is occupied by the BCS-type excitations. Both types of excitations have the width of 
the same order of magnitude transferring their energy and momentum to the normal quasiparticle 
excitations. We now estimate 7. For the entire Fermi level occupied by the normal state, the width 
is equal to 7 « iV(0) 3 T 2 /e{T) 2 , with the density of states N(0) ~ M*(T) ~ 1/T (see Eq. (2.12)). 
The dielectric constant e(T) ~ N(0), and the width 7 becomes 7 ~ T 



In our case, however, 



only a part of the Fermi surface within 9 C is occupied by the normal excitations. Therefore, the 
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number of states allowed for quasiparticles and for quasiholes is proportional to 6 C , and the factor 
T 2 is replaced by T 2 # 2 . Taking this into account, we obtain 7 ~ # 2 T ~ T(T — T c )/T c ~ (T — T c ). 
Here we have omitted the small contribution coming from the BCS-type excitations. That is why 
the width 7 vanishes at T = T c . Thus, the presented above analysis shows that in the pseudogap 
state at T > T c , the superconducting gap smoothly transforms into the pseudogap. The excitations 
of that area of the Fermi surface that has the gap are of the same width 7 ~ (T — T c ). The region 
occupied by the pseudogap is shrinking with increasing temperature. It is worth noting that the 
normal state resistivity p{T) behaves as p{T) oc T due to 7 ~ (T — T c ). Obviously at T > T* the 
behavior p(T) oc T remains valid up to temperatures T ~ Tf, and Tf can be as high as the Fermi 
energy, provided that FC occupies noticeable part of the Fermi volume, see Eq. (2.11). 

The temperature T no & e is determined mainly by the repulsive interaction being part of the pair- 
ing interaction A V A (pi,p 2 ). In its turn, the repulsive interaction can depend on the properties of 
materials such as the composition, doping, etc. Since the superconductivity is destroyed at T c , the 



ratio 2Ai/T c can vary in a wide range and strong 
follows from considerations given above j^" 



y depends upon the material's properties, as it 
j. The ratio 2Ai/T c can reach very high values. 



For instance, in the case of Bi 2 Sr 2 CaCu2Q6+<5 where the superconductivity and the pseudogap are 
considered to be of the common origin, 2Ai/T c is about 28, while the ratio 2A X /T* ~ 4, which is 
in agreement with the experimental data for various cuprates Note that Eq. (3.16) gives also 
good description of the maximum gap Ai in the case of the d-wave superconductivity, because the 
different regions with the maximum absolute value of Ai and the maximal density of states can be 

n 

considered as disconnected 60J. Therefore the gap in this region is formed by the attractive phonon 
interaction, which is approximately independent of the momenta. We can also conclude, that in the 
case of the s-wave pairing the pseudogap phenomenon is absent because there is no the sufficiently 
strong repulsive interaction. Thus, the transition from superconducting gap to pseudogap can take 
place only in the case of the d-wave pairing, so that the superconductivity is destroyed at T c ~ T no d e , 
with the superconducting gap being smoothly transformed into the pseudogap, which closes at some 
temperature T* > T c (3, Km Q • In the case of the s-wave pairing, we expect the absence of the 
pseudogap phenomenon in accordance with the observations (see e.g. j3] and references therein). 
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3.5. Dependence of T c on the doping 



We now turn to consideration of the maximum value of the superconducting gap Ai as a function 
of the number density x of the mobile charge carriers which is proportional to the doping. Using Eq. 
(2.15), we can rewrite Eq. (3.10) as follow 

(3.17) 

e F x FC 

Here we take into account that the Fermi level ep oc p F , the density x ~ p|./(2M£), and thus, ep oc x. 
We can reliably assume that T c oc Ai because the empirically obtained simple bell-shaped curve of 
T c (x) in the high temperature superconductors Q] should have only a smooth dependence upon x. 
Then, T c (x) in accordance with the data has the same bell-shaped form 6s| 

T c (x) oc P(xfc — x)x. (3.18) 

As an example of the implementation of the previous analysis, let us consider the main features 
of a room-temperature superconductor. The superconductor has to be a quasi two-dimensional 
structure like cuprates. From Eq. (3.10) it follows, that A x ~ 0ep oc f3/r 2 s . Noting that FCQPT in 
3D systems takes place at r s ~ 20 and in 2D systems at r s ~ 8 38], we can expect that Ai of 3D 
systems comprises 10% of the corresponding maximum value of 2D superconducting gap, reaching 



a value as high as 60 meV for underdoped crystals with T c = 70 69]. On the other hand, it is seen 
from Eq. (3.10), that Ai can be even large, Ai ~ 75 meV, and one can expect T c ~ 300 K in the 
case of the s-wave pairing as it follows from the simple relation 2T C ~ Ai. Indeed, we can safely take 
Sp ~ 300 meV, (3 ~ 0.5 and (pj —Pi)/p F ~ 0.5. Thus, a possible room-temperature superconductor 
has to be the s-wave superconductor in order to get rid of the pseudogap phenomena, which can 
tremendously reduce the transition temperature T c . The density x of the mobile charge carriers must 
satisfy the condition x < xpc and be adjustable to reach the optimal doping level x opt ~ xpc/2. 



3.6. The gap and specific heat near T c 

Now we turn to the calculations of the gap and the specific heat at the temperatures T — > T c . It 
is worth noting that this consideration is valid provided T* = T c . Otherwise the considered below 
discontinuity in the specific heat is smoothed out over the temperature range T* T c . For the sake 
of simplicity, we calculate the main contribution to the gap and the specific heat coming from FC. 
The function A X (T — ► T c ) is found from Eq. (3.14) by expanding the right hand side of the first 
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integral in powers of Ai and omitting the contribution from the second integral on the right hand 
side of Eq. (3.14). This procedure leads to the following equation 



Ax(T) ~ 3.4Tjl-^ . 



(3.19) 



Thus, the gap in the spectrum of the single-particle excitations has the usual behavior. To calculate 
the specific heat, the conventional expression for the entropy S [49J can be used 

dp 



S = -2|[/(p)ln/(p) + (l-/(p))ln(l-/(p))] 



(2tt) s 



where 



/(p) 



l + exp[£(p)/T] 
The specific heat C is determined by the equation 

Eo 



£(p)-^ + Af(T) 



(3.20) 



(3.21) 



dT 



+ 4^/ f(E)(l-f(E)) 

E 



4^ / f(E)(l -/(£)) 



J>2 







E 2 + TA 1 (T) 



dAt{T) 



dT 



d£ + 



E 2 + TA 1 (T) 



^Ax(T) 
dT 



(3.22) 



In deriving Eq. (3.22) we again used the variable £ and the densities of states N FC and N L , just 
as before in connection with Eq. (3.9), and employed the notation E - 



/£ 2 + Af(T). Eq. (3.22) 

predicts the discontinuity 5C = C s — C n in the specific heat C at T c because of the two last term in 
the square brackets on right hand side of Eq. (3.22). Here, C s and C n are the specific heat of the 
superconducting state and the normal one respectively. Using Eq. (3.19) to calculate the first term 
on the right hand side of Eq. (3.22), we obtain 

3 



5C(T C 



2tt 2 



(Pf-Pi)p F - 



(3.23) 



This is in contrast to the conventional result where the discontinuity is a linear function of T c . 5C(T C ) 
is independent of the critical temperature T c because as seen from Eq. (3.13) the density of states 
varies inversely with T c . Note that in deriving Eq. (3.23) we took into account the main contribution 
coming from FC. This term vanishes as soon as E —>■ and the second integral on the right hand 
side of Eq. (3.22) gives the conventional result. 

A few remarks are in order here. As we shall demonstrate in Section 9, Eq. (9.4), the specific 



heat of systems with FC behaves as CiT) oc JT/Tf. The specific heat discontinuity given by Eq. 
(3.23) is temperature independent. As a result, we obtain that 

5C(T C ) fT)( Pf - Pi ) 



C n (T c ) 



Pf 



(3.24) 
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In contrast to the conventional case of normal superconductors, when 5C(T c )/C n (T c ) = 1.43 |27j ]. 
it is seen from Eq. (3.24) that the ratio 6C(T c )/C n (T c ) is not the constant and can be very large 
provided that T c /Tf <C 1. 

4. THE DISPERSION AND LINESHAPE OF THE SINGLE-PARTICLE SPECTRA 

The newly discovered additional energy scale manifests itself as a break in the quasiparticle disper- 
sion near (50 — 70) meV, which results in a drastic change of the quasiparticle velocity 4fi| . 
Such a behavior is qualitatively different from what one could expect in a normal Fermi liquid. 
Moreover, this behavior can hardly be understood in the frames of either the Marginal Fermi Liquid 
theory or the quantum protectorate since there are no additional energy scales, or parameters, in 



these theories 
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HQ- 



One could suggest that this observed strong self-energy effect, leading to 



the new energy scale, is due to the electron coupling with collective excitations. But in that case one 



41] 



has to give up the quantum protectorate idea, which would contradict observations 0, 

As we have seen in Sections 2 and 3, Eqs. (2.12) and (3.7), the system with FC is characterized 
by two effective masses: M FC that is related to the single particle spectrum at lower energy scale, 
and Ml describing the spectrum at higher energy scale. These two effective masses manifest itself as 
a break in the quasiparticle dispersion, which can be approximated by two straight lines intersecting 
at the energy Eq, Eqs. (2.14) and (3.8). This break takes place at temperatures T ^ Tf when the 
system both in its superconducting state and normal one. This beahvior is in good agreement with 
the experimental findings |44j. It is pertinent to note that at T < T c , the effective mass M FC does 
not depend on the momenta pp , pf and Pi as it follows from Eqs. (3.7) and (3.10), 

2tt 

M*p C ~ (4.1) 
Thus, it is seen from Eq. (4.1) that M FC does not depend on x. This result is in good agreement 



with experimental facts 



(4.1) t 
flfl 



741 ] . The same is true for the dependence of the Fermi velocity 
vf = Pf/M fc on x, because the Fermi momentum pp ~ ^fn only slightly depends on the electron 
density n = n (l — x) [z3, tz^j • Here n is the single-particle electron density at the half-filling. 

Since Ao is the coupling constant defining the pairing interaction, for example phonon-electron 
interaction, it could be reasonable to expect that the break in the quasiparticle dispersion comes from 
the phonon-electron interaction. The phonon scenario could explain the persistence of the break at 
T > T c , since phonons are T-independent. On the other hand, it was shown that the quasiparticle 
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dispersion tends to recover to the conventional one-electron dispersion when the energy is well above 
the typical phonon energies The experimental observations do not show that the recovery to 



the one-electron dispersion takes place 44 1. 



The lineshape function L(q,cj) of the single-particle spectrum is a function of two variables. 
Measurements carried out at a fixed binding energy u = Uq, with uoq being the energy of a single- 
particle excitation, determine the lineshape L(q, to = ujq) as a function of the momentum q. We have 
shown above that M FC is finite and constant at T < T c . Therefore, at excitation energies to < Eq, 
the system behaves like an ordinary superconducting Fermi liquid with the effective mass given by 
Eq. (3.7) QHQ. At T c < T the low energy effective mass M FC is finite and is given by Eq. 
(2.12). Once again, at the energies u> < Eq, the system behaves as a Fermi liquid, the single- par ticle 



spectrum is well defined while the width of single-particle excitations is of the order of T 
This behavior was observed in experiments measuring the lineshape at a fixed energy 



20J2 




The lineshape can also be determined as a function of a;, L(q = q , u), at a fixed q = q . At small u, 
the lineshape resembles the considered above, and L(q = qo,oo) has the characteristic maximum and 
width. At energies u > Eq, the contribution coming form quasiparticles with the mass M£ become 
important, leading to the increase of L(q = qo,u). As a result, the function L(q = qo,u>) possesses 
the known peak-dip- hump structure [l?! directly defined by the existence of the two effective masses 
M FC and M£ j^J 30, We can conclude that in contrast to the Landau quasiparticles, these 
quasiparticles have a more complicated lineshape. 

To develop deeper quantitative and analytical insight into the problem, we use the Kramers- 
Kronig transformation to construct the imaginary part Im£(p, e) of the sin gle- particle self-energy 
S(p, e) starting with the real one ReS(p, e), which defines the effective mass jZS| 

1 (\_ 1 gReS \ / / gReS \ 

Here M is the bare mass, while the relevant momenta p and energies e obey the following strong 
inequalities: \p — Pf\/pf 1, and e/ep 1. We take ReS(p, e) in the simplest form which accounts 
for the change of the effective mass at the energy scale Eq: 



v,/ X M Fr ( Eq\ M Fr — Mr 

Re E(p,e) = -e— ^+(s ' ' 



E \ , B f Eq 



2 J V 2 



(4.3) 



M \ 2 J M 

Here 6(e) is the step function. Note that in order to ensure a smooth transition from the single- 
particle spectrum characterized by M FC to the spectrum defined by M£ the step function is to be 
substituted by some smooth function. Upon inserting Eq. (4.3) into Eq. (4.2) we can check that 

27 



inside the interval (—Eq/2, Eq/2) the effective mass can be estimated as M* ~ M FC , and outside the 
interval it is M* ~ M* L . By applyin g th e Kramers- Kronig transformation to ReE(p, e), we obtain 
the imaginary part of the self-energy |5jJ 

Im E(p.e) ~ e 2 — ^ + — ^ ^ e n 0/ , 



+ ^ln 
2 



e 2 -£ 2 /4 



^o74 



(4.4) 



We see from Eq. (4.4) that at e/E <C 1 the imaginary part is proportional to e 2 , at 2e/E ~ 1 
ImS ~ e, and at E /e <C 1 the main contribution to the imaginary part is approximately constant. 
This is the behavior that gives rise to the known peak-dip-hump structure. It is seen from Eq. (4.4) 
that when E — > the second term on the right hand side tends to zero and the single-particle 
excitations become better defined, resembling the situation in a normal Fermi-liquid, and the peak- 
dip-hump structure eventually vanishes. On the other hand, the so called renormaliztaion constant, 
or the quasiparticle renormalization factor, a(p) is given by [3] 

1 _ ,.8^) 



a(p) de 

At T < T c , as seen from Eqs. (4.3) and (4.5), the quasiparticle amplitude at the Fermi surface rises 
as the energy scale E decreases. It follows from Eqs. (3.8) and (3.18) that E ~ (xpc ~ x)/xpc- 
At T > T c , it is seen from Eqs. (4.3) and (4.5) that the quasiparticle amplitude rises as the effective 
mass Mp C decreases. As seen from Eqs. (2.12) and (2.15), Mp C ~ (pf — Pi)/pF ~ (xfc — x)/xfc- 
As a result, we can conclude that the amplitude rises as the level of doping increases, while the 
peak-dip-hump structure vanishes and the single-particle excitations become better defined in highly 
overdoped samples. At x > xpc-, the energy scale -E^ = and the quasiparticles are normal excitations 
of Landau Fermi liquid. It is worth noting that such a behavior was observed experimentally in highly 
overdoped Bi2212 where the gap size is about 10 meV [7^. Such a small size of the gap verifies that 
the region occupied by the FC is small since Eq/2 ~ Ai. Then, recent experimental data have shown 
;hat the Landau Fermi liquid does exist in heavily overdoped non-superconducting Lai.7Sr .3CuO4 
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8l|. 



5. FIELD-INDUCED LFL IN HEAVY ELECTRON LIQUID WITH FC 

In this Section we consider the behavior of the heavy-electron liquid with FC in magnetic fields, 
assuming that the coupling constant is nonzero Ao 7^ but infinitely small. As we have seen in 
Section 3, at T = the superconducting order parameter ft(p) is finite in the FC range, while 
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the maximum value of the superconducting gap Ai oc Ao is infinitely small. Therefore, any small 
magnetic field B ^ can be considered as a critical field and will destroy the coherence of n(p) and 
thus FC itself. To define the type of FC rearrangement, simple energy arguments are sufficient. On 
one hand, while the field is zero in the sample until then the state with FC is not destroyed. The 
energy gain AEb due to removing the FC state is AEb oc B 2 and tends to zero with B — > 0. On the 
other hand, occupying the finite range (pf — pi) in the momentum space, the function n (p) given by 
Eq. (2.8), or by Eq. (3.6), leads to a finite gain in the ground state energy 0]. Thus, a new ground 
state replacing FC should have almost the same energy as the former one. Such a state is given by 
the multiconnected Fermi spheres resembling an onion, where the smooth quasiparticle distribution 
function n (p) in the (pf — pi) range is replaced by a multiconnected distribution z/(p) 



KP) =J2 (P~ V2k-l)0{p2k - p). 



(5.1) 



k=\ 



Here the parameters Pi < p\ < p 2 < ■ ■ ■ < P2n < Pf are adjusted to obey the normalization condition: 



P2k+3 dp 
V2k-1 



P2fe+3 dp 

n {p) 

P2k-1 



(2tt)3 J P2k ^ u ^(2tt)3* 
The corresponding multiconnected distribution is shown in Fig. 4. 



v(x) 



X 2 X 3 X 4 



X 2n-l X 2n 



FIG. 4: The function f (p) for the multiconnected distribution replacing the function no(p) in the range 
(pf — pi) occupied by FC so that pi < pf < Pf- 

For definiteness, let us consider the most interesting case of a 3D system, while the consideration 
of a 2D system also goes along the same line. We note that the idea of multiconnected Fermi spheres, 
with production of new, interior segments of the Fermi surface, has been considered some time ago 



84L |85| . Let us assume that the thickness 5p of each interior block is approximately the same 
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Sp ~ P2k+i — P2k and Sp is denned by B. Using the Simpson's rule, we obtain that the minimum 
loss in the ground state energy due to formation of the blocks is about (5p) 4 . This result can be 
understood by considering that the continuous FC function no(p) delivers the minimum value to 
the energy functional i£[7i(p)], while the approximation of ^(p) by the steps of size 5p produces the 
minimum error of the order of {Sp) 4 . Taking into account that the gain due to the magnetic field is 
proportional to B 2 and equating both of the contribution we obtain 

5p oc 4B. (5.2) 

Thus, at T — > when B — > 0, the thickness 5p goes to zero as well, 5p — ► 0, while the NFL behavior 
of the FC state is replaced by the behavior of LFL with the Fermi momentum pf. It follows from Eq. 
(3.6) that pj > pf, while the number density x of itinerant electrons remains constant. As we shall 
see, this observation plays important role when considering the Hall coefficient Rh{B) as a function 
of B at low temperatures in the HF metals with FC. 

To calculate the effective mass M*(B) as function of the applied magnetic field B, we observe 
that at T = the application of finite magnetic field B splits the FC state into the Landau levels 
and suppresses the superconducting order parameter K (p) thus destroying the FC state. Therefore 



the LFL behavior is expected to be restored [8|, |86j. The Landau levels at the Fermi surface can 
be approximated by a single block whose thickness in momentum space is Sp. Approximating the 
dispersion of quasiparticles within this block by e(p) ~ (p — p F + bp) (p — pp)/M, we obtain that the 
effective mass M*(B) ~ M/{5p/pp). The energy loss AEpc due to rearrangement of the FC state 
related to this block can be estimated using the Landau formula (^tJ 

*E FO = /(£(p)-^(p)^. (5.3) 

The region occupied by the variation Sn(p) has the thickness Sp, while (e(p)— p) ~ (p—pf)pf/M*(B). 
As a result, we have AEpc ~ Sp 2 /M*(B). On the other hand, there is a gain AEb ~ 
(B 2 fig) 2 M* (B)pf due to the application of the magnetic field and coming from the Zeeman splitting. 
Equating AE B to AE FC and taking into account that in this case M*(B) oc 1/Sp, we obtain the 
following relation 

Sr) 2 1 

« 777~77^t? ex B 2 M*(B). (5.4) 



M*(B) (M*(B)) 3 
It follows from Eq. (5.4) that the effective mass M*(B) diverges as 



M\B) ex 7=—= - (5.5) 

V-D — -D c0 

30 



Here, B c q is the critical magnetic field which drives both a HF metal to its magnetic field tuned QCP 
and the corresponding Neel temperature toward T = ^]. We note that in some cases B c0 = 0, for 
example, the HF metal CeRu2Si2 shows neither evidence of the magnetic ordering, superconductivity 
down to the lowest temperatures nor the LFL behavior |8j|. Equation (5.5) shows that by applying 
the magnetic field B > B c0 the system can be driven back into LFL with the effective mass M*(B) 
depending on the magnetic field. This means that the following dependences are valid: for the 
coefficient A(B) oc (M*(B)) 2 Q, for the specific heat, C/T = ^ (B) oc M*(B), and for the magnetic 
susceptibility Xo(B) oc M*(B). The coefficient A(B) determines the temperature dependent part of 
the resistivity, p{T) = p + Ap, with p being the residual resistivity and Ap = A(B)T 2 . Since the 
coefficient is directly determined by the effective mass, we obtain from Eq. (5.5) that 

MB) oc jr^- (5.6) 

It is seen that the well-known empirical Kadowaki- Woods (KW) ratio 9], K = Aj^ ~ const, is 
fulfilled. At this point, we stress that the value of K may be dependent on the degeneracy number of 
quasiparticles as it was recently shown. As a result, the grand-KW- relation produces good description 
of the data for whole the ranee of debate HF s^Q. In tfce shiest case when the heavy 
electron liquid is formed by quasiparticles with spin 1/2 and the degeneracy number is 2, K turns 
out to be close to the empirical value j^, called as the KW relation [jj. Therefore, we come to the 
conclusion that by applying magnetic fields the system is driven back into the LFL state where the 
constancy of the Kadowaki- Woods ratio is obeyed. 

At finite temperatures, the system remains in the LFL state, but there exists a temperature 
T*(B), at which the NFL behavior is restored. To calculate the function T*(B) , we observe that 
the effective mass M* characterizing the single particle spectrum cannot be changed at T*(B) since 
there are no any phase transitions. In other words, at the crossover point, we have to compare the 
effective mass M*(T) defined by T*(B), Eq. (2.12), and that M*(B) defined by the magnetic field 
B, Eq. (5.5), M*(T) ~ M*(B) 

1 1 



M*(T) M*(B) 

As a result, we obtain 



ocT*(£)oc— — oc JB-B c0 . (5.7) 



T*(B)^^B~B c0 . (5.8) 

At temperatures T > T*(B), the system comes back to the NFL behavior with M* defined by Eq. 
(2.12), and the LFL behavior disappear. We can conclude that Eq. (5.8) determines the line in the 
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B — T phase diagram that separates the region of the B dependent effective mass from the region 
of the T dependent effective mass. At the temperature T*(B), there occurs a crossover from the T 2 
dependence of the resistivity to the T dependence. It follows from Eq. (5.8), that the heavy electron 
liquid at some temperature T can be driven back into the Landau Fermi-liquid by applying a strong 
enough magnetic field (B — B c0 ) oc (T*(£>)) 2 . We can also conclude, that at finite temperature 
T < T*(B), the heavy electron liquid shows a more pronounced metallic behavior at the elevated 
magnetic field B since the effective mass decreases (see Eq. (5.5)). The same behavior of the effective 
mass can be observed in the Shubnikov — de Haas oscillation measurements. We conclude that one 
obtains a unique possibility to control the essence of the strongly correlated liquid by magnetic fields 
which induce the change of the NFL behavior to the LFL liquid behavior. 

Let us briefly consider the case when the system is very near FCQPT being on the ordered side 
and therefore Sppc — (Pf ~Pi)/Pf <C 1- Since 5p oc M*(B), it follows from Eqs. (5.2) and (5.5) that 



where a c is a constant, which is expected to be of the order of a unit, a c ~ 1. At bigger magnetic field 
B, the value of bpjpp becomes comparable with 5pFC) an d the distribution function ^(p) vanishes 
being replaced by the conventional Zeeman splitting. As a result, we are dealing with the heavy 
electron liquid located on the disordered side of FCQPT. As we shall see in Section 9, the behavior 
of this system is quite different from that of the system with FC. It follows from Eq. (5.9) that 
relatively weak magnetic field B cr 



removes the system from the ordered side of the phase transition provided that 8pfc *C I. 

6. MAGNETIC-FIELD INDUCED LANDAU FERMI LIQUID IN HIGH-T C METALS 

The LFL theory has revealed that the low-energy elementary excitations of a Fermi liquid look 
like the spectrum of an ideal Fermi gas. These excitations are described in terms of quasiparticles 
with an effective mass M*, charge e and spin 1/2. The quasiparticles define the major part of 
the low-temperature properties of Fermi liquids. As we have shown in Section 5, at temperatures 
T < T*(B), the LFL behavior of the heavy electron liquid with FC is recovered by the application of 
magnetic field B larger than the critical field B c suppressing the superconductivity. Thus, the heavy 
electron liquid with FC can be viewed as LFL induced by the magnetic field. In such a state, the 




(5.9) 



(5pfc) 2 




(5.10) 
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Wiedemann-Franz (WF) law and the Korringa law are held and the elementary excitations are LFL 
quasiparticles. Our consideration is valid for relatively weak magnetic fields when the contributions 
coming from the magnetic field are proportional B 2 and the Zeeman splitting is much smaller then 
the Fermi momentum. 

It was reported recently that in the normal state obtained by applying a magnetic field greater 
han the upper critical filed B c , in a hole-doped cuprates at overdoped concentration (Tl 2 Ba 2 Cu0 6+( ,) 
~y\ and at optimal doping concentration (Bi 2 Sr 2 Cu06+<5) js^]; there are no any sizable violations 
of the WF law. In the electron- doped copper oxide superconductor Pr a91 LaCeo.o9Cu0 4 _j / (T c =24 
K) when superconductivity is eliminated by a magnetic field, it was found that the som-lattice 
relaxation rate 1/Ti follows the T{F = constant relation, known as the Korringa law |2J|, down 
to temperature of T = 0, K Q. At higher temperatures and appfied magnetic fields of 15.3 T 
perpendicular to the Cu0 2 plane, 1/T\T as a function of T is a constant below T = 55 K. At 300 
K > T > 50 K, \jT{T decreases with growing T [9j|. Recent measurements for strongly overdoped 
non-superconducting La i 7S m sCuOd have shown that the resistivity p exhibits T 2 behavior, and the 
WF law holds perfectly |£Qj, . Since the validity of the WF and the Korringa laws are the robust 
signature of LFL, these experimental facts demonstrate that the observed elementary excitations 
cannot be distinguished from the Landau quasiparticles. This imposes strong constraints for models 
describing the hole-doped and electron-doped high-temperature superconductors. For example, in 



the cases of a Luttinger liquid 



93j . spin-charge separation (see e.g. (3)) an d i n some solutions of 



t — J model |9_5j a violation of the WF law was predicted. 

As any phase transition, FCQPT is related to the order parameter, which induces a broken 
symmetry. It was shown in Section 3 that the relevant order parameter is the superconducting order 
parameter n(p), which is suppressed by the critical magnetic field B c , when B 2 C ~ Af. If the coupling 
constant Ao — > 0, the critical magnetic field B c — ■> will destroy the state with FC converting the 
strongly correlated Fermi liquid into LFL. The magnetic field plays the role of the control parameter 
determining the effective mass M*(B) as it follows from Eq. (5.5). 

If Ao is finite, the critical field is also finite, and Eq. (5.5) is valid at B > B c . In that case the 
system is driven back to LFL and has the LFL behavior induced by the magnetic field. Then, the 
low energy elementary excitations are characterized by M*(B) and cannot be distinguished from 
Landau quas iparticles. As a result, at T — > 0, the WF law is held in accordance with experimental 
facts [ill IqoI]. On the hand, in contrast to the LFL theory, the effective mass M*(B) depends on the 
magnetic field. 
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Equation (5.5) shows that by applying a magnetic field B > B c the system can be driven back into 
LFL with the effective mass M*(B) which is finite and independent of the temperature. This means 
that the low temperature properties depend on the effective mass in accordance with the LFL theory. 
In particular, the resistivity p(T) as a function of the temperature behaves as p(T) = po + Ap(T) 
with Ap(T) = AT 2 , and the factor A oc (M*(B)) 2 . Taking into account that in the case of the 
high-T c superconductors B c0 is expected to be zero, we obtain from Eq. (5.5) that 

7oV^B = const. (6.1) 

Here 70 = C/T with C is the specific heat. Taking into account Eqs. (5.6) and Eq. (6.1), we obtain 

70 ~ A(B)VB. (6.2) 

At finite temperatures, the system remains LFL, but there is the crossover from the LFL behavior 
to the non- Fermi liquid behavior at temperature T*(B) oc y/~B. At T > T*(B), the effective mass 
starts to depend on the temperature M* oc 1/T, and the 



resistivity possesses the non-Fermi liquid 



behavior with a substantial linear term, Ap(T) oc T |29l I30L l83j . Such a behavior of the resistivity 
was observed in the cuprate superconductor Tl 2 Ba 2 Cu0 6+( 5 (T c < 15 K) [96| . At B=10 T, Ap(T) is a 
linear function of the temperature between 120 mK and 1.2 K, whereas at B = 18 T, the temperature 
dependence of the resistivity is consistent with Ap(T) = AT 2 in the same temperature range j^. 
In LFL, the nuclear spin-lattice relaxation rate 1/Ti is determined by the quasiparticles near the 



Fermi level whose population is proportional to M*T, so that 1/T±T oc M* is a constant |9l|, l92 |. 
When the superconducting state is removed by the application of a magnetic field, the underlying 
ground state can be seen as the field induced LFL with effective mass depending on the magnetic 
field. As a result, the rate 1/T"i follows the T{T = constant relation, that is the Korringa law is 
held. Unlike the behavior of LFL, as it follows from Eq. (5.5), 1/T{T oc M*(B) decreases with 
increasing the magnetic field at T < T*(B). At T > T*(B), we observe that \jT\T is a decreasing 
function of the temperature, \jT{T oc M* oc 1/T. These observations are in a good agreement with 
the experimental facts 9^. Since T*(B) is an increasing function of the magnetic field, see Eq. 
(5.8), the Korringa law retains its validity to higher temperatures at elevated magnetic fields. We 
conclude, that at temperature T < T*(B ) and bigger magnetic fields B > B the system shows 
a more pronounced metallic behavior, since the effective mass decreases with increasing B (see Eq. 
(5.5)). Such a behavior of the effective mass can be observed in the de-Haas van Alphen-Shubnikov 
studies, \jT{T and the resistivity measurements. These experiments can shed light on the physics 
of high-T c metals and reveal relationships between high-T c metals and heavy-electron metals |97| |. 
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The existence of FCQPT can also be revealed experimentally because at densities x > Xpc, or 
beyond the FCQPT point, the system should be LFL at sufficiently low temperatures [86]. Recent 
experimental data have shown that this liquid exists in the heavily overdoped non-superconducting 
compound Lai. 7 Sr .3 CuC-4 HQ- 

It is remarkable that up to T = 55 K the resistivity exhibits the 
T 2 behavior with no additional linear term, and the WF law is verified to within the experimental 



resolution 



While at elevated temperatures, a strong deviations from the LFL behavior are 



observed. We anticipate that in this case the system can be again driven back to the LFL behavior 
by the application of sufficiently strong magnetic fields. 

Thus, the mentioned above striking measurements, which were used in studies of the nature of 
the high-T c superconductivity, suggest that FCQPT and the emergence of the novel quasiparticles 
with effective mass strongly depending on the magnetic field and temperature and resembling the 
Landau quasiparticles are qualities intrinsic to the electronic system of the high-T c superconductors. 



7. APPEARANCE OF FCQPT IN DIFFERENT FERMI LIQUIDS 

It is widely believed that unusual properties of correlated liquids observed in the high-temperature 
superconductors, heavy- fermion metals, 2D 3 He and etc., are determined by quantum phase tran- 
sitions. Therefore, immediate experimental studies of relevant quantum phase transitions and of 
their quantum critical points are of crucial importance for understanding the physics of the high-T c 
superconductivity and HF systems. In case of the high-T c superconductors, these studies are diffi- 
cult to carry out, because all the corresponding area is occupied by the superconductivity. On the 
other hand, recent experimental data on different highly correlated Fermi liquids, when the system 
in question is approaching FCQPT from the disordered side, can help to illuminate both the nature 
of this point and the control parameter, by which this point is driven. We shall call Fermi liquids 
approaching FCQPT from the disordered side as highly correlated ones to distinguish them from 
strongly correlated liquids which have undergone FCQPT. Detailed explanations on this point are 
given in Section 8. 

Expe— facts on Mgh-densi* 2D 3 H e HQ show that the effective ma ss diverges wee. the 
density, at which 2D 3 He liquid begins to solidify, is approached 2lJ. Then, a sharp increase of the 
effective mass in a metallic 2D electron system was observed, when the density tends to the critical 

n 

density of the metal-insulator transition point, which occurs at sufficiently low densities [341] . Note, 
that there is no ferromagnetic instability in both Fermi systems and the relevant Landau amplitude 
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F a > -1 



34, 



99J], in accordance with the almost localized fermion model [10 



Now let us consider the divergence of the effective mass in 2D and 3D highly correlated Fermi 
liquids at T = 0, when the density x approaches FCQPT from the side of normal Landau Fermi 
liquid, that is from the disordered phase. First, we calculate the divergence of M* as a function of 
the difference (x — xpc) m case of 2D Fermi liquid. For this purpose we use the equation for M* 



obtained in 



381 ] . where the divergence of the effective mass M* due to the onset of FC in different 



Fermi liquids (such as 2D and 3D electron and 3 He liquids) was predicted. At x — > xfc, the effective 
mass M* can be approximated as 

1^1 1 j J v(q(y)) ydydg 

W~M + IW J[i- R(q(y),u = 0, g) X o(q(y), w = 0)] 2 v / W r l '" J 



Here we adopt the notation ppy2(l — y) = q(y) with q(y) being the transferred momentum, M is 
the bare mass, u is the frequency, v(q) is the bare interaction, and the integral is taken over the 
coupling constant g from zero to its real value go- In Eq. (7.1), both xo(<7,^) and R(q,u) are the 
linear response function of a noninteracting Fermi liquid and the effective interaction respectively. 
They define the linear response function of the system in question 

X{q,u,g) = - — r — r. (7.2) 

1 - R{q,uj,g)xo{q,uj) 

In the vicinity of the charge density wave instability, occurring at the density x^d w , the singular part 
of the function x 1 on the disordered side is of the well-known form (see e.g. 2j) 

X^iq^iQ) - a(x cdw -x) + b(q - q c f + c(g - g), (7.3) 

where a, b and c are constants and q c ~ 2pp is the wavenumber of the charge density wave order. 
Upon substituting Eq. (7.3) into Eq. (7.1) and integrating, the equation for the effective mass M* 
can be presented in the following form ju[ |2| 

J- = l °— (74) 

M* M ^/x - x cdw ' v ' ' 

with C being some positive constant. The behavior of the effective mass as a function of the electron 
number density i in a silicon MOSFET is shown in Fig. 5. The fitting parameters are x c d w = 
0.7 x 10- n cm- 2 , C = 2.14 x 10~ 6 cm" 1 and x FC = 0.9 x lO^cm" 2 3- lt is seen from Fi S- 5 that 
Eq. (7.4) describes rather good the data. 

It is seen from Eq. (7.4) that M* diverges at some point x fc_ referred to as the critical point, at 
which FCQPT occurs, as a function of the difference (x — xfc) 



HQ 



M B 

— ~A + , (7.5) 

M x-x FC 
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FIG. 5: The ratio M* /M in a silicon MOSFET versus the carrier number density x. The filled squares 
denote the Shubnikov - de Haas oscillations experiment al da ta, the data obtained by the application of a 
parallel magnetic field are shown by the filled circles 

where A and B are constants. It follows from the derivation of Eqs. (7.4) and (7.5) that their 
forms are independent of the bare interaction v(q), which effects however A, B and xpc values. This 
result is in agreement with Eq. (2.7) which exhibits the same type of divergence independent of the 
interaction. Therefore both of these equations are also applicable to 2D 3 He liquid or to another 
Fermi liquid. It is also seen from Eqs. (7.4) and (7.5) that FCQPT precedes the formation of 
charge- density waves. As a consequence of this, the effective mass diverges at high densities in case 
of 2D 3 He, and at low densities in case of 2D electron systems, in accordance with experimental facts 
34, 9^. Note, that in both cases the difference (x — xpc) has to be positive, because x approaches 
Xpc when the system is on the disordered side of FCQPT with the effective mass M*(x) > 0. Thus, 
in considering the 2D 3 He liquid we have to replace (x — xpc) by (xpc ~ x ) on the right hand side 
of Eq. (7.5). In case of a 3D system, at x — > xfc, the effective mass is given by 

1 1 p F } J v{q{y))ydydg 



M* M + A^ll [l-R(q( y ),u; = 0,g) Xo (q(y),u; = 0)f ^ 

The comparison of Eq. (7.1) and Eq. (7.6) shows that there is no fundamental difference between 
these equations, and along the same way we again arrive at Eqs. (7.4) and (7.5). The only difference 
between 2D electron systems and 3D ones is that in the latter FCQPT occurs at densities which are 
well below those corresponding to 2D systems. For bulk 3 He, FCQPT cannot probably take place 
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FIG. 6: The ratio M*/M in a 2D 3 He versus the fluid number density x inferred from the heat capacity 
measurements and the magnetization measurements 99j. The solid line represents M* /M = A + ^ Xp ^_ x ^ , 
with A=1.09, B=1.68 nm~ 2 and xpc = 5.11 nm -2 . 

since it is absorbed by the first order solidification 99|. The apparent divergence of the effective 
mass M*(x) obtained in measurements on 2D 3 He [99] is shown in Fig. 6. This is in good agreement 
with the divergence given by Eqs. (7.5) and (2.7). 



8. HEAVY FERMION METALS WITH HIGHLY CORRELATED ELECTRON LIQUID 

In HF metals with strong electron correlations, quantum phase transitions at zero temperature 
may strongly influence the measurable quantities up to relatively high temperatures. These quan- 
tum phase transitions have recently attracted much attention because the behavior of HF metals is 
expected to follow universal patterns defined by the quantum mechanical nature of the fluctuations 
taking place at quantum critical points (see e.g. 0,Q). Only recently, there appeared experimen- 
tal facts which deliver experimental grounds to understand the nature of quantum phase transition 
producing the universal behavior of HF metals. It has been demonstrated that at low temperatures 
the main properties of HF metals such as the magnetoresistance, resistivity, specific heat, magneti- 
zation, susceptibility, volume thermal expansion, etc, strongly depend on temperature T and applied 
magnetic field B. As a result, these properties can be controlled by placing these metals at the 
special point of the field-temperature B — T phase diagram. In the LFL theory, considered as the 
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main instrument when investigating quantum many electron physics, the effective mass M* of quasi- 
particle excitations determining the thermodynamic properties of electronic systems is practically 
independent of temperature and applied magnetic fields. Therefore, the observed anomalous behav- 
ior is uncommon and can hardly be understood within the framework of the conventional LFL theory 
based on the notion of quasiparticles. As a result, it is necessary to use theories that are based on 
the Landau concept of the order parameter which is introduced to classify phases of the state of 
matter. These theories connect the anomalous behavior with critical fluctuations of the magnetic 
order parameter. These fluctuations suppressing the quasiparticles are attributed to CQPT taking 
place when the system in question approaches its QCP. As it was noted in Introduction, the universal 
behavior is only observable if the electron system of HF metal is very near QCP, for example, when 
the correlation length is much larger than the microscopic length scales. In the case of CQPT, the 
physics is dominated by thermal and quantum fluctuations of the critical state, which is characterized 
by the absence of quasiparticles. It is believed that the absence of quasiparticle-like excitations is 
the main cause of the NFL behavior and other types of the critical behavior in the quantum criti- 
cal region. However, theories based on CQPT fail to explain the experimental observations related 
to the divergence of the effective mass M* at the magnetic field tuned QCP, the specific behavior 
of the spin susce ptibi l ity and its scaling properties, the therma l expansion behavior, etc, see e.g. 



14. 



87 



102, 



103, 



104, 



105 



106 



107, 



108, 



109 



110, 



The LFL theory rests on the notion of quasiparticles which represent elementary excitations 
of a Fermi liquid. Therefore these are appropriate excitations to describe the low temperature 
thermodynamic properties. The inability of the LFL theory to explain the experimental observations 
which point to the dependence of M* upon the temperature T and applied magnetic field B may 
lead to the conclusion that the quasiparticles do not survive near QCP, and one might be further led 
to the co nclu s ion t hat the heavy electron does not retain its integrity as a quasiparticle excitation 



(see e.g. 



102, 



10 



ml 



HQ). 



The mentioned above inability to explain the behavior of HF metals at QCP within the framework 



of theories based on CQPT may also lead to the con e 



of the order parameter fails as well, see e.g. 
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usion that the other important Landau concept 



107 



112 



113j |. Thus, we are left without the 



most fundamental principles of many body quantum physics while a great deal of interesting NFL 
phenomena related to the anomalous behavior and the experimental facts collected in measurements 
on the HF metals remain out of reasonable theoretical explanations. 

On the other hand, it is the very nature of HF metals that suggests that their unusual properties 
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are defined by a quantum phase transition related to the unlimited growth of the effective mass at 
its QCP. Moreover, a di vergence to infinity of the effective electron mass was observed at a magnetic 
field- induced QCP In Section 2, we have demonstrated that such a quantum phase 

transition is to be FCQPT, an essential feature of which is the divergence of the effective mass M* 
at its QCP, see Eq. (2.7). 

8.1. Highly correlated heavy electron liquid 

When a Fermi system approaches FCQPT from the disordered phase it remains the Landau Fermi 
liquid with the effective mass M* strongly depending on the distance r = (x — £fc), temperature 
T and a magnetic field B. This state of the system, with M* essentially depending on T, r and B, 
resembles the strongly correlated liquid described in Section 2. However, in contrast to the strongly 
correlated liquid, there is no energy scale Eq given by Eq. (2.14) and the system under consideration 
is the Landau Fermi liquid at sufficiently low temperatures with the effective mass M* oc 1/r (see Eqs. 
(2.7) and (7.5)). Therefore this liquid can be called the highly correlated liquid which is obviously 



to have uncommon properties |2(1 |3l|. Again, we use the heavy electron liquid model to study the 
universal behavior of the HF metals at low temperatures. As it was mentioned in Section 2, it is 
possible, since we consider processes related to the power-low divergence of the effective mass. This 
divergence is determined by small momenta transferred as compared to momenta of the order of the 
reciprocal lattice cell, and the contribution coming due to the lattice structure can be ignored. Thus, 
we may usefully ignore the complications related to the lattice and get rid of the specific peculiarities 
of a HF metal regarding the medium as homogeneous heavy electron isotropic liquid. 

The effective mass M* of quasiparticle excitations controlling the density of states determines 
the thermodynamic properties of electronic systems. To study the behavior of the effective mass 
M*(T, B) as a function of temperature and magnetic field, we use the Landau equation determining 
M*(T, B). In the case of a homogeneous liquid, at finite temperatures and low magnetic fields, this 
equations reads [27 1 

1 1 , [ PfPi p / ,dn ai (p 1 ,T,B) dp 1 



2^ J — 3-^Wpf,pi 



M*(T,B) M p F ^^'^ dp 1 (2tt) 3 ' 

Here F ayCri (pp, pi) is the Landau amplitude depending on the momenta p and spins a, pf is the Fermi 
momentum, M is the bare mass of an electron and n a (p, T) is the quasiparticle distribution function. 
Since HF metals are predominantly three dimensional (3D) structures we treat the homogeneous 
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heavy electron liquid as a 3D liquid also. For the sake of simplicity, we omit the spin dependence 
of the effective mass since in the case of a homogeneous liquid and weak magnetic fields, M*(T,B) 
does not noticeably depend on the spins. The quasiparticle distribution function is of the form 



n a (p,T) = \ 1 + exp 



(e(p,T) - fi a ) 
T 



-i 



(8.2) 



where e(p, T) is determined by Eq. (2.2). In our case, the single-particle spectrum does not noticeably 
depend on the spin, while the chemical potential may have a dependence due to the Zeeman splitting. 
We will show explicitly the spin dependence of a physical value when this dependence is of importance 
for understanding. 

Replacing n CT (p, T, B) by n a (p, T, B) = Sn a (p, T, B) + n a (p, T — 0, B — 0) where Sn a (p, T, B) = 
n CT (p, T, B) - n a (p, T = 0, B = 0), Eq. (1) takes the form 

M M M rp FPl ^ d5n ai ( Pl ,T,B) dp, 

+ 3"Z^ /— ^ CT ,<xi(PF,Pij tt 77^- (.8.6) 



M*(T,B) M*(x)p 2 F ^J p F w*'^' 9 Pl (2tt)3- 

We assume that the heavy electron liquid is near FCQPT, therefore the distance r is small so that 
M/M*(x) <C 1, as it is seen from Eq. (2.7). In the case of normal metals with the effective mass 
of the order of a few bare electron masses, M/M*(x) ~ 1, and up to temperatures T ~ 100 K, the 
second term on the right hand side of Eq. (8.3) is of the order of T 2 / ji 2 and is much smaller than the 
first term. Thus, the system in question demonstrates the LFL behavior with the effective mass being 
practically independent of temperature, that is the corrections are proportional to T 2 . One can check 
that the same is true when magnetic field up to B ~ 30 T is applied. Near the critical point xfc, 
when M/M*(x — > x F c) — > 0, the behavior of the effective mass changes drastically, because the first 
term on the right hand side of Eq. (8.3) vanishes and the second term determines the effective mass 
itself rather than small corrections to M*(x) related to T and B. In that case, Eq. (8.3) becomes a 
homogeneous equation and determines the effective mass as a function of B and T. As we will see, 
Eq. (8.3) describes both the NFL behavior and the LFL one with the presence of quasiparticles. In 
contrast to the conventional Landau quasiparticles these are characterized by the effective mass that 
strongly depends on both the magnetic field and the temperature. 

Let us turn to a qualitative analysis of solutions of Eq. (8.3) when x ~ xfc- We start with the 
case when T = and B is finite. The application of magnetic field leads to the Zeeman splitting of 
the Fermi surface and the difference Sp between the Fermi surfaces with "spin up" and "spin down" 
becomes 5p = p F — pp ~ fi B BM*(B)/p F with fi B being the Bohr magneton. Upon taking this into 
account, we observe that the second term in Eq. (8.3) is proportional to (5p) 2 oc (fi B BM*(B)/p F ) 2 , 

41 



and Eq. (8.3) takes the form |sl l8fiUll4 



M M , \y B BM\B))\ (g 4) 



M*(B) M*(x) p% 
where c is a constant. Note that the effective mass M*(B) depends on x as well and this dependence 
disappears at x — Xpc- At the point x = xfcj , the ter m M/M*(x) vanishes, Eq. (8.4) becomes 



homogeneous and can be solved analytically |3l|, |86J, 1114 1 




Here B c0 is the critical magnetic field which drives both a HF metal to its magnetic field tuned QCP 
and the corresponding Neel temperature toward T = 8]. We recall that in some cases B c0 = 0, for 
example, the HF metal CeRu2Si2 shows neither evidence of the magnetic ordering, superconductivity 
down to the lowest temperatures nor the LFL behavior ^&\ . 

Equation (8.5) shows the universal power low behavior of the effective mass which does not depend 



on t 



18, 



re inter-particle interaction. We illustrate this behavior by calculations using a model functional 



IS | 

2 



with the inter-particle interaction 

= go expR3o| P |) ( 

IpI 

We normalized the effective mass by M, M* = M*(B)/M, temperature To and magnetic field H 
by the Fermi energy e° F , T = T /e° F , H = (hbB)/e Fi and use the dimensionless coupling constant 
g = (goM)/(2/pi 2 ) and (3 = /3qPf- FCQPT takes place when the parameters reach their critical 
values, (3 = b c and g = g c . On the other hand, FCQPT takes place when effective mass M* —* oo. 
This condition allows to relate b c and g c [1$, 

0(1 + 6 c )exp(-6 c )[6 c cosh(6 c ) - sinh(6 c )] = 1. 

c 

It follows from this equation that the critical point of FCQPT can be reached by changing go if /3q 
and pp are fixed, or changing pp if /3 and go are fixed, etc. For simplicity, we shall change g to reach 
FCQPT and investigate the properties of the system behind the critical point. 

In Fig. 7 we present our calculations (triangles and squares) of the magnetic-field dependence 
of the effective mass at the critical point of FCQPT. At (3 = b c = 3, FCQPT takes place when 
g = g c = 6.7167. It is seen from Fig. 7 that the calculated power-low divergence of the effective mass 
is in accordance with Eq. (8.5). 
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FIG. 7: Calculated magnetic-field dependent effective mass M* at fixed temperatures. Arrows indicate 
the temperatures. The vertical axis represents the normalized effective mass M*. The horizontal axis is 
the normalized magnetic field H. Details of normalization are explained in the text. Solid lines represents 
M*(H) oc H- 2 / 3 . 

At densities x > xpc, M*(x) is finite and we are dealing with the conventional Landau quasipar- 
ticles provided that the magnetic field is weak, so that M*(x)/M*(B) 1 with M*(B) given by Eq. 
(8.5). In that case, the second term on the right hand side of Eq. (8.4) is proportional to (BM*(x)) 2 
and represents small corrections. In the opposite case, when M*(x)/M*(B) ^> 1, the heavy electron 
liquid behaves as at QCP. Since in the LFL regime the main thermodynamic properties of the system 
is determined by the effective mass, it follows from Eq. (8.5) that we obtain a unique possibility to 
control the magnetoresistance, resistivity, specific heat, magnetization, volume thermal expansion, 
etc. At this point, we note that the large effective mass leads to the high density of states provoking 
a large number of states and phase transitions to emerge and compete with one another. Here we 
assume that these can be suppressed by the application of a magnetic field and concentrate on the 
thermodynamical properties. 

To consider the qualitative behavior of M*(T) at elevated temperatures, we simplify Eq. (8.3) 
by omitting the variable B and simulating the influence of the applied magnetic field by the finite 
effective mass entering the denominator of the first term on the right hand side of Eq. (8.3). This 
effective mass becomes a function of the distance r, M*(r), which is determined by both B and 
(x — xfc)- If the magnetic field vanishes the distance is r = (x — xfc)- We integrate the second term 
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over the angle variable, then over p% by parts and substitute the variable p\ by z, z = (e(pi) — fij/T. 
In the case of the flat and narrow band, we use the approximation (e(pi) — /i) — ~ Pf)/M*(T) 

and finally obtain 

MM/" 00 1 r 1 /" 1 

, : /rTnN = , - , : + a / F(Pf, + az)) dz — a F(pf,Vf({ — otz)) dz. (8.8) 

M*(T) M*(r) Jo "l + e 2 Jo ;; l + e* v ; 

Here the notations are used: F ~ Md{F 1 p 2 ) / 'dp, the factor a = TM*(T)/p 2 F = TM*(T)/(T k M*(r)), 

T k = p 2 F /M*(r) The Fermi momentum is defined as e(pp) = \i. We first assume that a < 1. Then 

omitting terms of the order of exp(— 1/a), we expand the upper limit of the second integral on the 

right hand side of Eq. (8.8) to oo and observe that the sum of the second and third terms represents 



an even function of a. These are the typical expressions with 



ermi-Dirac functions as integrands. 



They can be calculated using standard procedures (see e.g. [115I ]). Since we need only an estimation 
of the integrals, we represent Eq. (8.8) as 

M M (TM*(T)\ 2 /TM*(T)\ 4 , , 

+ «iUrni4 +••• (8-9) 



M*(T) M*(r) \T k M*(r)J \T k M*(r 
Here a\ and 02 are constants of the order of units. Equation (8.9) can be considered as a typical 
equation of the LFL theory with the only exception being the effective mass M*(r) which strongly 
depends on the distance r = x — xfc > ® and diverges at r — > 0. Nonetheless, it follows from Eq. 
(8.9) that when T — > the corrections to M*(r) start with the T 2 terms provided that 

M fTM*(T)\ 2 T 2 

» Urr^y (8-10) 



M*{r) \T k M*(r) ) T 2 ' 
and the system exhibits the LFL behavior. It is seen from Eq. (10) that when r — ► 0, M*(r) — > 00, 
and the LFL behavior disappears. The free term on the right hand side of Eq. (8.8) vanishes, 
M/M*(r) — > 0, and Eq. (8.8) in itself becoming homogeneous determines the value and universal 
behavior of the effective mass. 

At some temperature T\ T k , the value of the sum on the right hand side of Eq. (8.9) is 
determined by the second term. Then Eq. (8.10) is not valid, and upon ke eping only the second 
term in Eq. (8.9) this can be used to determine M*(T) in a transition region 114L lllfil ] 

M\T) oc JL (8.11) 

The variation as T~ 2 / 3 exponent with the temperature growth deserves a comment. Equation (8.11) 
is valid if the second term in Eq. (8.9) is much larger than the first one, that is 

T 2 M . 

n » wFy (ai2) 
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and this term is grater than the third one, 

T M*(r) 

T t «jFw-r L (8 ' 13) 

Obviously, both Eqs. (8.12) and (8.13) can be simultaneously satisfied if M/M*(r) <C 1 and T is 
finite. The range of temperatures, over which Eq. (8.11) is valid, shrinks to zero as soon as r — > 
because Tk — > 0. Thus, if the system is very near QCP, x — > xfc, it is possible to observe the behavior 
of the effective mass given by Eq. (8.11) in a wide range of temperatures provided that the effective 
mass M* (r) is diminished by the application of high magnetic fields, that is, the distance r becomes 
larger under the action of B. When r is finite the T~ 2 / 3 behavior can be observed at relatively high 
temperatures. To estimate the transition temperature Ti(B), we observe that the effective mass is a 
continuous function of the temperature, thus M*(B) ~ M*(Ti). Taking into account Eqs. (8.7) and 
(8.11), we obtain T^B) oc B. 




FIG. 8: Calculated entropy S(T) as a function of temperature at fixed magnetic field H = 0.01. Solid line 
represents S(T) at the transition region T = T\(B) (left panel). At low temperatures T < 0.007, the system 
exhibits the LFL behavior, S(T) oc T. The other solid (right panel) line represents S(T) oc T 1 / 2 . 

Then, at elevated temperatures, the system enters into a different regime. The coefficient a 
becomes a ~ 1, the upper limit of the second integral in Eq. (8.8) cannot be expanded to oo, and 
odd terms come into play. As a result, Eq. (8.9) is no longer valid, but the sum of both the first 
integral and the second one on the right hand side of Eq. (8.8) is proportional to M*(T)T. Upon 
omitting the first term M/M*(r) and approximating the sum of the integrals by M*(T)T, we solve 
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Eq. (8.8) and obtain 

M*(T) oc J=. (8.14) 



We illustrate the above consideration by numerical calculations (shown by filled squares in Fig. 
8) based on the model functional (8.6). In Fig. 8, we show the evolution of the low temperature 
entropy from the transition region T ~ T\(B) to S(T) oc \/T at higher temperatures upon applying 
magnetic field H = 0.01. The calculated behavior of S(T)/T oc M*(T) is in accord with Eq. (8.14). 
Details of the normalization T and H are given at Fig. 7. 

Thus, we can conclude that at higher temperatures when x ~ xfc the system exhibits three types 
of regimes: the LFL behavior at a < 1, when Eq. (8.10) is valid; the M*(T) oc T~ 2//3 behavior 
and S(T) oc M*(T)T oc T 1 / 3 , when Eqs. (8.12) and (8.13) are valid; and the 1/VT behavior 
of the effective mass at a ~ 1, while the entropy S(T) oc M*(T)T oc a/T, and the specific heat 
C{T) = T(dS(T)/dT) oc VT. 

In the absence of magnetic field, calculated evolutions of M*(T), S(T) and C(T) based on the 
model functional (8.6) are shown by filled squares in Fig. 9, Fig. 10, and Fig. 11 respectively. 




FIG. 9: Calculated effective mass M*(T) as a function of temperature. Solid line represents M*(T) oc 1/VT. 

Let us estimate the quasiparticles width 7(T). Within the framework of the LFL theory it is given 
by H 

7 ~ |r| 2 (M*) 3 T 2 , (8.15) 

where T is the particle-hole amplitude. In the case of a strongly correlated system with its large 
density of states related to the huge value of the effective mass, the amplitude T cannot be approxi- 
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FIG. 10: Calculated entropy S(T) as a function of temperature. Solid line represents S(T) oc yT 
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FIG. 11: Calculated specific heat C(T) as a function of temperature. Solid line represents C(T) oc yT 



mated by the bare particle interaction but can be estimated within the ladder approximation which 
gives |T| ~ l/(ppM*(T)) [20j. As a result, we obtain that in the LFL regime j(T) oc T 2 , in the 
rp-2/3 re gj me ry(T) oc T 4 / 3 , and in the 1/VT regime 7(T) oc T 3 / 2 . We observe that in all the cases 
the width is small as compared to the quasiparticle characteristic energy which is of the order of T, 
so the notion of a quasiparticle is meaningful. We can conclude that when the heavy electron liquid 
is near the QCP of FCQPT, being on the disordered side, its low energy excitations are quasiparticle 
excitations with the effective mass M*(B, T). At this point we note that at x — > xpc-, the quasiparti- 
cle renormalization factor a(p) remains finite and approximately constant, and the divergence of the 
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effective mass given by Eq. (2.7) is not related to vanishing a(p) 117j . Thus the notion of the quasi- 
particles is preserved and these are the relevant excitations when considering the thermodynamical 
properties of the heavy electron liquid. 



8.2. Resistivity of heavy fermion metals 

Since the resistivity, p(T) = po + Ap(T), is directly determined by the effective mass, because 
the coefficient A(B,T) oc (M*(B,T)) 2 j88|, the above mentioned temperature dependences can be 
observed in measurements of the resistivity of HF metals. 

At T <C Tx, the system in question demonstrates the LFL regime, the divergence of the effective 
mass at x — > xpc is described by Eq. (8.5) and the coefficient A(B) diverges as 

Thus, the resistivity, p(B) oc A(B)T 2 , as a function of magnetic field diverges as given by Eq. (8.16). 
As a function of temperature, the resistivity behaves as Api = CxT 2 / (B — 5 c0 ) 4//3 oc T 2 . The second is 
NFL regime which is determined by Eq. (8.11) and characterized by Ap 2 = c{T 2 j i^T 2 ^) 2 oc T 2 / 3 . At 
T > Tx(B), the third NFL regime is given by Eq. (8.14) and represented by Ap% = C3T 2 /(\/T) 2 oc T. 
Here cj, C2, C3 are constants. It is remarkable that all temperature dependences corr esponding to thes e 
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104. 
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109]. 



regimes were observed in measurements on the HF metals CeCoIn 5 and YbAgGe 
If we consider the ratio Ap 2 / Api oc ((B — S c o)/T) 4//3 , we come to a very interesting conclusion that 
the ratio is a function of only the variable (B — B c q)/T, thus representing the scaling behavior. This 
result is in 



8.3. Magnetic susceptibility 

The magnetic susceptibility is proportional to the effective mass, \ oc M*, with M* given by Eq. 
(8.5). Therefore, at T <^ Tx, 

X (B) oc M*(B) oc (B - £ c0 r 2/3 , (8.17) 
while the static magnetization Mb{B) is given by 

M B (B) oc BM*(B) oc(B- B c0 ) 1/3 . (8.18) 
At T ^> Tx, as it follows from Eq. (8.14), Eq. (8.17) has to be rewritten as 

x(T) oc M*(T) oc —f=. (8.19) 
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The observed behavior of x{B) and Mb(B) and the behavior of x(T) are in accord with results of 
measurements on CeRu 2 Si 2 with the critical field B c0 — > |stI ] . 

Consider the state of the system when r — > 0. Its properties are determined by magnetic fields B 
and temperature T because there are no other parameters to describe such a state. At the transition 
temperatures T ~ Ti(B), the effective mass depends on both T and B, while at T T\(B), the 
system is LFL with the effective mass being given by Eq. (8.5), and at T > Ti(B), the mass is 
defined by Eq. (8.11). Instead of solving Eq. (8.3), it is possible to construct a simple interpolation 
formula to describe the behavior of the effective mass over all regions, 

M *^' T ) = m r wu T( ytw - (8 ' 20) 
d{B - Bco) 2 ' 6 + c 2 f{y)T I l- i 

Here, f(y) is a universal monotonic function of y = (T/(B — B c q)) 2 ^ such that f(y ~ 1) = 1, and 

f(y < 1) = 0. It is seen from Eq. (8.20) that the behavior of the effective mass can be represented 

by a universal function of only one variable y if the temperature is measured in the units of the 

transition temperature Ti(B), and the effective mass is measured in the units of M*(B) given by Eq. 

(8.5). This representation describes the scaling behavior of the effective mass. As seen from Eqs. 

(8.18) and (8.20), the scaling behavior of the magnetization can be represented in the same way, 

provided the magnetization is normalized by the saturated value at each field given by Eq. (8.5) 

M b{B,T) 1 

M B (B) 0C l + c S /(y)y' ^ 

where C3 is a constant. It is seen from Eq. (8.21), that magnetization is a monotonic function of 

y. Upon using the definition of susceptibility, x — dM B /dB, we come to the conclusion that the 

susceptibility also exhibits the scaling behavior and can be presented as a universal function of only 

one variable y, provided it is normalized by the saturated value at each field given by Eq. (8.17) 

xCgiS g 1 ,o c J(v) + ydf(y)/dy 

It is of importance to note that the susceptibility is not a monotonic function of y because the 
derivative is the sum of two contributions with different behavior. The second contribution on 



the right hand side of Eq. (8.22) gives the susceptibility a maximum |8j, Ill4l Ill6 |. As shown in 



the susceptibility and magnetization |114j and with the facts observed in measurements 



Fig. 12, the above behaviors of the magnetization and susceptibility are in accord with numerical 
calculations o: 
on CeRu 2 Si 2 [87 1. 

It is seen from Fig. 12 that for finite B, the curve describing x{B, T) acquires a maximum at some 
temperature Tp. This behavior is in good ag ^twithe.pe^al^Ha.d — with 
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FIG. 12: Normalized magnetic susceptibility x(B,T)/x(B,Tp) (left panel) and normalized magnetization 
Mb(B,T)/Mb(B,Tp) (right panel) for CeRu2Si2 in magnetic fields 0.20 mT (squares), 0.39 mT (triangles), 



and 0.9 I mT (circles), plowed against normalized temperature 7'/7;» Jg^, w iere L)> is the Irnipemlure at. 
peak susceptibility. The solid curves trace the calculated universal behavior 11 J |. 



le LFL behavior 
1141 . Note that 



the suppression of the divergent NFL terms ~ T~ 2 ^ 3 in x(B,T) and recovery of t 
at static magnetic fields in which the Zeeman energy splitting ppB exceeds T 
magnetization Mb(B,T) does not exhibits a maximum. In Fig. 12, the temperature is normalized 
by Tp, the susceptibility is normalized by its peak height x{B,Tp) and the magnetization by the 
saturated value at each field. 



8.4. Magnetoresistance 

Using the results just presented, we consider the behavior of magnetoresistance (MR) 

Pmr(B,T)- _ , (8.23) 

as a function of magnetic field B and T. Here the resistivity, p(B, T) = p + Ap(B, T) + Ap mr (B), 
is measured at the magnetic field B and temperature T. We assume that the contribution Ap mr (B) 
coming from the magnetic field B can be treated within the low field approximation and given by 
the well-known Kohler's rule, 



Ap mr {B) (Bp(0,e D 



2 



*n-MwJ' (a24) 

with 0£) is the Debye temperature, B is the characteristic field and is a constant. Note, that 

the low field approximation implies that Ap mr (B) p(0,T) = p(T). We also assume that that 

50 



temperature is not too low so that po < Ap(B = 0, T), while B > 5 c o. Substituting Eq. (8.24) into 
Eq. (8.23), we find that 

, R ^ c(M*(E,T)) 2 T 2 + Ap mr .{B) - c(M*(0,T)) 2 T 2 

Pmr{B,T) . (8.25) 

Here M*(B,T) denotes the effective mass which now depends on both the magnetic field and the 
temperature, and c is a constant determining the temperature dependent part of the resistivity, 
c(M*(B,T)) 2 T 2 = Ap(B,T). 

Consider MR given by Eq. (8.25) as a function of B at some temperature T = To. At low magnetic 
fields when To > T\(B) oc B, the main contribution to MR comes from Ap mr (B) since the effective 
mass depends mainly on temperature. Therefore, the ratio \M*(B,T) — M*(0,T)\/M*(0,T) <C 1, 
the main contribution is given by Ap mr (B), and MR is an increasing function of B. When B 
becomes so large that T\(B) ~ To, the difference (M*(B,T) — M*(0,T)) becomes negative, and MR 
as the function of B reaches its maximum value at T\{B) ~ To. We recall that T\{B) determines 
the crossover from T 2 dependence of the resistivity to the T dependence. At elevated B when 
Tx(B) > T , the ratio (M*(B,T) - M*(0, T))/M*(0, T) ~ -1 and MR becomes negative being a 
decreasing function of B. 

Consider now MR as a function of T at some B . At T < Tx(B Q ), we have LFL. At low tem- 
peratures T < Ti(B ), it follows from Eqs. (8.5) and (8.14) that M*(B )/M*(T) < 1, and MR 
is determined by the resistivity p(0,T). Note, that Bq has to be comparatively high to ensure the 
inequality, M*(B )/M*(T) < 1. As a result, p mr (B ,T -> 0) ~ -1, because Ap mr (B)/p(0, T) < 1. 
Differentiating the function p mr (Bo, T) with respect to Bo we can check that its slope becomes steeper 
as Bq is decreased, being proportional oc (B —B c o)~ 7 / 3 . At T ~ Ti(Bq), MR possesses a node because 
at this point the effective mass M*(B ) ~ M*(T), and p(B , T) ~ p(0, T). We can conclude that the 
crossover from the T 2 resistivity to the T resistivity, which occurs at T ~ Ti(Bq), manifests itself 
in the transition from negative MR to positive MR. At T > Ti(Bq), the main contribution to MR 
comes from Ap mr (B ), and MR reaches its maximum value. Upon using Eqs. (8.14) and (8.25) and 
taking into account that at this point T oc (B — B c0 ), we obtain that the maximum value p^ r (-B ) 
of MR is p™ r (B ) oc 1/(B — B c0 ). Thus, the maximum value is a decreasing function of B . At 
T\(Bq) T, MR is a decreasing function of the temperature. At these temperatures, MR becomes 
small comparatively to its maximum value p™ r (B ) because \M*(B,T) - M*(0, T)\/M*(0, T) < 1 
and Ap mr (B )/p(T) <L 

The recent paper |l05| reports the measurements of the CeCoIn 5 resistivity in magnetic fields. 
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The both transitions from negative MR to positi ve M R with increasing T and from positive MR 
to negative one with increasing B were observed jl05| . Thus, the above observed behavior of MR 
is in good agreement with the experimental facts. We believe that an additional analysis of the 



data 



1051 ] can reveal that the crossover from T 2 dependence of the resistivity to the T dependence 



occurs at T oc (B — B c0 ) (B — B c0 )^ 3 , this analysis could reveal the above described supplementary 
peculiarities of MR as well 



9. HF METALS WITH STRONGLY CORRELATED ELECTRONIC LIQUID 

As we have seen in Section 2, at T = when r = (x—xpc) ~^ 0, the effective mass M*(r) — > oo and 
eventually beyond the critical point x = xpc the distance r becomes negative making the effective 
mass negative as follows from Eq. (2.7). As it was shown in Section 2, the system is to undergo 
FCQPT. Therefore behind the critical point xfc of this transition, the quasiparticle distribution 
function represented by the step function does not deliver the minimum to the Landau functional 
E[n(p)]. As a result, at x < xpc the quasiparticle distribution is determined by Eq. (2.8) to 
search the minimum of a functional, which determines the quasiparticle distribution function no(p) 
delivering the lowest possible value to the ground state energy E. It was shown in Section 3 that the 
relevant order parameter /t(p) = yVio(p)(l — ^o(p)) is a t the same time the order parameter of the 
superconducting state when the maximum value Ai of the superconducting gap is infinitely small. 
Thus, this state cannot exist at any finite temperatures and driven by the parameter x: at x > xpc 
the system is on the disordered side of FCQPT, while at x < xfc-, the system is on the ordered 
side. In Section 2 we have shown that this ordered state has a strong impact on the properties of 
the system at finite temperatures T <C Tf. the effective mass M*(T) diverges as 1/T, see Eq. (2.12), 
and the electronic system with FC is characterized by the energy scale Eq given by Eq. (2.14). As a 
result, we can consider FCQPT as the phase transition, which separates the highly correlated heavy 
electron liquid from the strongly correlated one. It was shown in Section 8, that the highly correlated 
liquid at T — > and x > xpc behaves as LFL, therefore, FCQPT separates the regions of LFL and 
strongly correlated liquids. Obviously, the strongly correlated electron liquid demonstrates the NFL 
behavior down to zero temperature. 

At < T < T/, the function no(p) determines the entropy Snfl{T) of the heavy electron 
liquid in its NFL state. Inserting into Eq. (2.3) the function no(p), one can check that behind the 
point of FCQPT there is a temperature independent contribution S (r) ~ (pf — pi)/pF ~ \ r \/ x FC, 
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where r = x — xfc- Another specific contribution is related to the spectrum e(p) which insures the 
connection between the dispersionless region (pf —pi) occupied by FC and the normal quasiparticles 
located at p < pi and at p > pf, and therefore it is of the form e(p) oc (p — pf) 2 ~ (pi — p) 2 . Such 
a form of the spectrum can be verified in exactly solvable models for systems with FC and leads to 
the contribution of this spectrum to the specific heat C ~ ' /Tf |l7[. Thus at < T <C Tf, the 
entropy can be approximated as 

S NFL (T) ~ S (r) + a^- + b^, (9.1) 

with a and b are constants. The third term on the right hand side of Eq. (9.1) comes from the 
contribution of the temperature independent part of the spectrum e(p) and gives a relatively small 
contribution to the entropy. 

The calculated evolution of S(T), C(T) and M*(T) based on the model functional (8.6) are shown 
by filled symbols in Fig. 13, Fig. 14, and Fig. 15 respectively. The calculations were carried out for 
g = 7, 8, 12 and /3 = b c = 3, while the critical value g = g c = 6.7167. 




FIG. 13: Calculated entropy S(T) as a function of temperature. Solid lines represent S(T) given by Eq. 
(9.1). 

It is seen from Fig. 13 that So(r) increases when the system being on the ordered side moves 
away from FCQPT. 

Obviously, the term S'o(r) does not contribute to the specific heat. As a result the specific heat 
demonstrates the anomalous behavior C(T) oc y/T as it is seen from Fig. 14. As to the effective 
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FIG. 14: Calculated specific heat C(T) as a function of temperature. Solid lines represent C(T) oc vT- 




— i 1 1 1 1 1 1 1 1 — 1 1 — 

0,00 0,02 0,04 0,06 0,08 0,10 



T 

FIG. 15: Calculated effective mass M*(T) as a function of temperature. Solid lines represent M*(T) oc 
ai/T + a 2 VT + a 3 . 

mass M*(T) oc S(T)/T, it demonstrates the divergence M*(T) oc 1/T in accordance with Eq. (2.5) 
and contains all terms defining the behavior of the entropy, as shown in Fig. 15. 

9.1. The Griineisen ratio and Hall coefficient in heavy fermion metals; T—B phase diagram 

The temperature independent term So(r) determines the specific NFL behavior of the system. For 
example, the existence of the temperature independent term So(r) of the entropy can be illuminated 



54 



by calculating the thermal expansion coefficient ct(T) [71. Ill8l |. which is given by 



= 1 fdQogV)\ x fdjS/a 



^ (9 2) 

3 V dT J p 3K \ dx j T K ' ' 

Here, P is the pressure and V is the volume. The compressibility K is not expected to be singular 
at FCQPT and in systems with FC, because FC is attached to the Fermi level, and it m oves along 



as [a(x) changes, while the compressibility K = dfi/d(Vx) is approximately constant 119l |. Inserting 
Eq. (9.1) into Eq. (9.2), we find that 

M* T 

a FC {T) ~ a ~ -g^. (9.3) 

Here, is a number independent of temperature. When deriving Eq. (9.3) we keep only the main 
contribution coming from Sq(t). On the other hand, the specific heat 

As a result, the Griineisen ratio T(T) diverges as 

At this point, we consider how the behavior of the effective mass given by Eqs. (2.12) and 
(8.14) correspond to experimental observations. It was recently observed that the thermal expansion 
coefficient a (T) / T measured on CeNi2 Ge2 shows a 1/ VT divergence over two orders of magnitude in 
the temperature range from 6 K down to at least 50 mK, while measurements on YbRh2(Sio.95Geo.os)2 
demonstrate that a/T oc 1/T [6j], contrary to the LFL theory which yields a(T)/T cx M* ~ const. 
Since the effective mass depends on T, we obtain that the 1/VT behavior, Eq. (8.14), is in excellent 
agre ement with the result for the former system |7j, and the 1/T behavior, Eq. (2.12), predicted in 
[ll8j corresponds to the latter HF metal. 

We see that at < T Tf, the heavy electron liquid with FC behaves as if it were placed 
at QCP. In fact it is placed at the quantum critical line x < xfc, that is the critical behavior is 
observed at T — > for all x < xpc- At T — > 0, the heavy electron liquid undergoes a first-order 
quantum phase transition because the entropy is not a continuous function: at finite temperatures 
the entropy is given by Eq. (9.1), while S(T = 0) = 0. Therefore, the entropy undergoes a sudden 
jump 5S = So(r) in the zero temperature limit. We make up a conclusion that due to the first 
order phase transition, the critical fluctuations are suppressed at the quantum critical line and the 
corresponding divergences, for example the divergence of T(T), are determined by the quasiparticles 
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rather than by the critical fluctuations as one could expect in the case of CQPT, see e.g. [4]. Note 
that according to the well known inequality, 5Q < T5S, the heat 5Q of the transition from the 
ordered phase to the disordered one is equal to zero, because 5Q < So(r)T — > at T — > 0. 

To study the B — T phase diagram of the heavy electron liquid with FC, we consider the case when 
the NFL behavior arises by the suppression of the antiferromagnetic (AF) phase upon applying a 
magnetic field B, for example, as it takes place in the HF metals YbRh 2 Si2 and YbRh 2 (Sio.95Geo.os)2 
The AF phase is represented by the heavy electron LFL, with the entropy vanishing as T — > 0. 
For magnetic fields exceeding the critical value B c0 at which the Neel temperature T^{B — ► B c0 ) — > 
the weakly ordered AF phase transforms into weakly polarized heavy electron LFL. As it was 
discussed in Section 5, at T = the application of the magnetic field B splits the FC state occupying 
the region (pf — Pi) into the Landau levels and suppresses the superconducting order parameter 
«(p) destroying the FC state. Such a state is given by the multiconnected Fermi sphere, where the 
smooth quasiparticle distribution function %(p) in the {pj—Pi) range is replaced by a multiconnected 
distribution z/(p), see Fig. 4. Therefore the LFL behavior is restored being represented by the weakly 
polarized heavy electron LFL and characterized by quasiparticles with the effective mass M*(B) given 



by Eq. (5.5). At elevated temperatures T > T*(B — B c0 ) oc y/B — B c0 , the NFL state is restored 
and the entropy of the heavy electron liquid is given by Eq. (9.1). This behavior is displayed in the 
T — B phase diagram shown in Fig. 16. 

In accordance with experimental facts we assume that at relatively high temperatures T/Tjyo ~ 1 
the AF phase transition is of the second order 40], where T^o is the Neel temperature in the 
absence of the magnetic field. In that case, the entropy and the other thermodynamic functions are 
continuous functions at the transition temperature Tn(B). This means that the entropy of the AF 
phase Saf(T) coincides with the entropy of the NFL state given by Eq. (9.1), 

S AF (T - T N (B)) = S NFL (T - T N (B)). (9.6) 

Since the AF phase demonstrates the LFL behavior, that is Saf{T — > 0) — > 0, Eq. (10) cannot be 
satisfied at sufficiently low temperatures T < T crit due to the temperature-independent term S (r), 
see Eq. (9.1). Thus, the second order AF phase transition becomes the first order one at T = T crit 
as it is shown in Fig. 16. At T = 0, the critical field B c q, at which the AF phase becomes the heavy 
LFL, is determined by the condition that the ground state energy of the AF phase coincides with the 
ground state energy _E[n (p)] of the heavy LFL, that is the ground state of the AF phase becomes 
degenerated at B = B c0 . Therefore, the Neel temperature T^{B — > B c0 ) — > 0, and the behavior of 
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FIG. 16: T — B phase diagram of the heavy electron liquid. The T^(B /B c q) curve represents the field 
dependence of the Neel temperature. Line separating the antiferromagnetic (AF) and the non-Fermi liquid 
(NFL) state is a guide to the eye. The black dot at T = T cr a marked by the arrow is the critical temperature, 
at which the second order AF phase transition becomes the first one. At T < T cr n, the thick solid line 
represents the field dependence of the Neel temperature when the AF phase transition is of the first order. 
The NFL state is characterized by the entropy Snfl given by Eq. (9.1). The line separating the NFL state 
and the weakly polarized heavy electron Landau Fermi Liquid (LFL) is T*(B/B c q — 1) oc ^/B/B c q — 1. 

the effective mass M*(B > S c0 ) is given by Eq. (5.5), that is M*(B) diverges when B — > B cQ . 

We note that the corresponding quantum and thermal critical fluctuations vanish at T < T crit 
because we are dealing with the first order AF phase transition. We can also reliably conclude 
that the critical behavior observed at T — > and B — > B c q is determined by the corresponding 
quasiparticles rather than by the critical fluctuations accompanying second order phase transitions. 
When r — > the heavy electron liquid approaches FCQPT from the ordered phase. Obviously, 
Tcrit — > at the point r = 0, and we are led to the conclusion that the Neel temperature vanishes at 
the point when the AF second order phase transition becomes the first order one. As a result, one 
can expect that the contributions coming from the corresponding critical fluct uatio ns can only lead 
to the logarithmic corrections to the Landau theory of the phase transitions 115j . and the power 



low critical behavior is again defined by the corresponding quasiparticles. Thus, we conclude that 
the Landau paradigm based on the notions of the quasiparticles and order parameter is applicable 
when considering the heavy electron liquid. 

Now we are in position to consider the recently observed jump in the Hall coefficient at B — > B c0 
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in the zero temperature limit 120|. At T = 0, the application of the critical magnetic field B c q 
suppressing the AF phase (with the Fermi momentum p^p — Pf) restores the LFL with the Fermi 
momentum Pf > Pf- At B < B c q, the ground state energy of the AF phase is lower then that of the 
heavy LFL, while at B > B c q, we are dealing with the opposite case, and the heavy LFL wins the 
competition. At B — B c0 , both AF and LFL have the same ground state energy being degenerated. 
Thus, at T = and B = B c0 , the infinitesimal change in the magnetic field B leads to the finite 
jump in the Fermi momentum because the distribution function becomes multiconnected, see Fig. 4, 
while the overall Fermi volume remains constant. That is, the number of itinerant electrons does not 
change. In response the Hall coefficient Rh{B) oc 1/x oc l/p 3 f undergoes the corresponding sudden 
jump. Here we have assumed that the low temperature Rh(B) can be con sidered as a measure of 



the Fermi volume and, therefore, as a measure of the Fermi momentum 



12(|. As a result, we obtain 



^ = ^°''j - 1 + 3 P -l^ K 1 + d^. (9.7) 
R H (B = B cQ + 5) p F x fc 

Here 5 is an infinitesimal magnetic field, Sq{t)/xfc is the entropy per one heavy electron, and d is a 

constant, d ~ 5. It follows from Eq. (9.1) that the abrupt change in the Hall coefficient tend s to z ero 
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when r — > and vanishes when the system in question is on the disordered side of FCQPT 

Now consider the magnetic susceptibility which is proportional to the effective mass given by Eq. 
(5.5). Therefore, at T T*(B), the magnetic susceptibility given by Eq. (5.8) is of the form 

X (B)<xM*(B)<x 7 =L==, (9.8) 

while the static magnetization M(B) is given by jq| 



M(B) (x^fB- B c0 . (9.9) 
As seen from Fig. 17, the field dependence of M(B) given b y Eq . (9.9) is in good agreement with 



the data obtained in measurements on YbRh 2 (Sio.95Ge .o5)2 [1Q2[. We can also conclude that Eqs. 
(8.21) and (8.22) determining the scaling behavior of the effective mass, static magnetization and the 
susceptibility are also valid in the case of strongly correlated liquid, but the variable y is now given by 



y = T/y/B — B c0 , while the function f(y) can be dependent on (pf — pi)/pF- This dependence comes 
from Eq. (2.12). As a result, we can obtain that at T < T*(B), the factor dp/dT oc A(B)T behaves 
as A(B)T oc T/(B — B c q), and at T > T*(B), it behaves as A(B)T oc 1/T. These o bserv ations are 



in good agreement with the data obtained in measurements on YbRh 2 (Si .95Ge .o5)2 |102 1 . 

We note that, as in the case of the highly correlated liquid, the susceptibility x(5,T) of the 
strongly correlated liquid is not a monotonic function of y and possesses a maximum as a function 
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FIG. 17: Magnetization M{B) shown by filled squares as a function of magnetic field B 1021 ] . The curve 
represents the field dependence of M(B) = clmVB given by Eq. (9.9) with ayi is a costnat. 

of the temperature because the derivative dM(B)/dB is the sum of two contributions. As it was 
shown in Section 5, the well-known empirical Kadowaki- Woods ratio jjj, K = A/75 2 — const , is also 
obe yed in the case of the strongly correlated liquid. These results are in good agreement with facts 
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As an application of the above consideration we study the T — B phase diagram for the HF 
metal YbRli2Si2 shown in Fig. 18. The LFL behavior is characterized by the effective mass 
M*(B) which diverges as 1 / \JB — B c0 [10]. We can conclude that Eq. (5.5) gives good description 
of this experimental fact, and M*(B) diverges at the point B — > B c0 with T^(B = B c0 ) = 0. It 
is seen from Fig. 18, that the line separating the LFL state and NFL can be approximated by the 



function cV B — B c0 with c being a parameter. Taking into account that the behavior of YbRh 2 Si2 
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103, 122], we can conclude that in 



strongly resembles the behavior of YbRh 2 (Sio.95Geo.o5)2 
the NFL state the thermal expansion coefficient a(T) does not depend on T and the Griineisen ratio 
as a function of temperature T diverges p|. We are led to the conclusion that the entropy of the 
NFL state is given by Eq. (9.1). Taking into account that at relatively high temperatures the AF 
phase transition is of the second order [10], we predict that at lower temperatures this becomes the 
first order phase transition. Then, the described behavior of the Hall coefficient Rh{B) is in good 
agreement with experimental facts 120 ]. 



Thus, we can conclude that the T — B phase diagram of the heavy electron liquid with FC is in 
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FIG. 18: T — B phase diagram for YbRh 2 Si2 JlOj, L102] . The T/v curve represents the field dependence of the 



Neel temperature. Line separating the antiferromagnetic (AF) and the non-Fermi liquid (NFL) state is a 
guide to the eye. The NFL state is characterized by the entropy Snfl given by Eq. (9.1). Line separating 



the NFL state and LFL is T*(B — B c q) = c\/B — B c q, with c being an adjustable factor. 

good agreement with the experimental T — B phase diagram obtained in measurements on the HF 
metals YbRh 2 Si2 and YbRb^Sio.gsGeo.os^- 



9.2. Heavy fermion metals very near FCQPT on the ordered side 

Let us consider the case when 5ppc — ipf ~ Pi)/Pf 1 and the electronic system of HF metal 
is very near FCQPT being on the ordered side. As we have seen in Section 5, Eq. (5.10), the 
application of magnetic field (B — B c0 )/B c0 > B c removes the system from the ordered side placing 
it on the disordered one. As a result at T < Ti(B), the effective mass M*(B) is given by Eqs. 
(8.5) and (8.11), while the corresponding resistivity is described in Section 8.2. In the absence of 
magnetic field or at Tf ^> T > Ti(B), the system demonstrates the NFL behavior, the effective 
mass M*(T) is given by Eq. (2.12), and the entropy is given by Eq. (9.1). While the magnetic 
susceptibility x(T) oc M*(T) oc l/T, the thermal expansion coefficient ot(T) is T-independent being 
determined by Eq. (9.3) and making the Griineisen ratio divergent, as it follows from Eq. (9.5). 
Note, that the specific heat behaves as C(T) oc y/T in both cases: when the electronic system is on 
the ordered side or on the disordered side of FCQPT, see Fig. 11 and Fig. 15. It follows from Eq. 
(2.12) that 7(T) cx T. Therefore, the temperature dependent part Ap(T) of the resistivity behaves 
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as Ap(T) oc j(T) oc T. Thus, the system demonstrates the NFL regime, Ap(T) oc T, when it is 
either in the highly correlated or in the strongly correlated regimes. 

At some temperature T c , the system can undergo the superconducting phase transition. In contrast 
to the conventional superconductors where the discontinuity 5C(T C ) in the specific heat at T c is a 
linear function of T c , SC(T C ) is independent of the critical temperature T c . As seen from Eqs. (3.23) 
and (3.24), both the disconti nuity 5C(T C ) and the ratio SC(T C ) /C n (T c ) can be very large as compared 



to the conventional case 54 
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Recent experiments show that the electronic system of the HF metal CeCoIn 5 can be considered as 
the system located near FCQPT and containing FC. Indeed, under the application of magnetic field 
it be haves as the highly correlated liquid (see Section 8) with the effective mass A(B) oc (B — S c o) _4//3 
jl05^ . as it is seen from Eq. (8.16). While in the NFL regime, ct(T) oc const and the Griineisen ratio 



diverges 



Ui, - Eq , (9.3) and (9.5), respectively. Estpat.ons of SpFC based on the evasion 



ie o btained value of Sp^c allows 
12311 observed at T r = 2.3 K in 



of the magnetic susceptibility show that 5ppc — 0.044 [123J. T 
to explain the relatively bi g va lue of the discontinuity 5C(T C 
measurements on CeCoIns |l2fij . Thus, we can conclude that B cr ~ 0.003, as it follows from Eq. 
(5.10), and the HF metal CeCoIn 5 can be considered as placed near FCQPT. 

10. DISSYMMETRICAL TUNNELLING IN HF METALS AND HIGH-T C SUPERCON- 
DUCTORS 

Experiments on the HF metals explore mainly their thermodynamic properties. It would be 
desirable to probe the other properties of the heavy electron liquid such as the probabilities of 
quasiparticle occupations, whi ch are not directly linked to the density of states or to the behavior 



of the effective mass M* [l26j|. Scanning tunnelling microscopy (STM) being sensitive to both the 
density of states and the probabilities of quasiparticle occupations is an ideal technique for studying 
such effects at quantum level. 

The tunnelling current / through the point contact between two ordinary metals is proportional to 
the driving voltage V and to the squared modulus of the qua ntum mechanical transition amplitude 
t multiplied by the difference A r i(0)A" 2 (0)(ni(p, T) — n^ip., T)) \l2*\ . Here n(p,T) is the quasiparticle 
distribution function and N(0) is the density of states of the corresponding metal. On the other 
hand, the wave function calculated in the WKB approximation and defining t is proportional to 
(A" 1 (0)A r 2 (0)) -1 / 2 . As a result, the density of states is dropped out and the tunnelling current is 
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independent of A^i(0)A^(0). Upon taking into account that at T — ► the distribution n(p, T —>■ 0) — > 
rip(p), where rip(p) is the step function 6(p — pp), one can check that within the LFL theory the 
differential tunnelling conductivity cr<i(V) — dl/dV is a symmetric function of the voltage V. 

In fact, the symmetry of &d(V) holds provided that so called particle-hole symmetry is preserved 
as it is within the LFL theory. Therefore, the existence of the UdiY) symmetry is quite obvious and 
common in the case of metal-to-metal contacts when these metals are ordinary metals and in their 
normal or superconducting states. 

Now we turn to a consideration of the tunnelling current at low temperatures which in the case 
of ordinary metals is given by |l 2t| 



I(V) = 2\t\ 2 J [n F (z -fi)- n F (z - fi + V)\ dz. (10.1) 

We use an atomic system of units: e = m = h = 1, where e and m are electron charge and mass, 
respectively, and the energy z belongs to the interval Eq given by Eq. (2.14) 

fi - 2T < z < fi + 2T. (10.2) 

Since temperatures are low, we approximate the distribution function of ordinary metal by the step 
function rip. It follows from Eq. (10.1) that quasiparticles with the energy z, fi — V < z < fi, 
contribute to the current, while &d{V) — 2\t\ 2 is a symmetrical fu nctio n of V. In the case of the 
heavy electron liquid with FC, the tunnelling current is of the form jl2fi| | 

I{V) = 2 J [n (z - fi)- n F (z - fi + V)] dz. (10.3) 

Here we have replaced the distribution function of ordinary metal by Uq that is the solution of Eq. 
(2.8). We have also normalized the transition amplitude \t\ 2 such that \t\ 2 = 1. Assume that V 
satisfies the condition, \V\ < 2T, while the current flows from the HF metal to the ordinary one. 
Quasiparticles of the energy z, fi — V < z, contribute to I(V), and the differential conductivity is 
givem by the relation er^V) ~ 2n (z ~ fi — V). If the sign of the voltage is changed, the direction 
of the current is also changed. In that case, quasiparticles of the energy z, fi + V > z, contribute to 
/(V), and the differential conductivity <J d (—V) ~ 2(1 — n (z ~ p, + V)). The dissymmetrical part 
Acd(V) = (crd(—V) — °~d(V)) of the differential conductivity is of the form 

Aa d {V) ~ 2[1 - {n {z -fi~V)- n (z - fi ~ -V))]. (10.4) 

It is worth noting that according to Eq. (10.4) we have A<Jd{V) = if the considered HF metal is 
replaced by an ordinary metal. Indeed, the effective mass is finite at T — > 0, then n (T —>■ 0) —>■ rip 
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being given by Eq. (2.4), and 1 — n(z — \i ~ V) = n(z — fi ~ — V). One might say that the 
dissymmetrical part vanishes due to the particle-hole symmetry. On the other hand, there are no 
reasons to expect that (1 — n (z — \i ~ V) — n (z — // ~ —V)) = 0. Thus, we conclude that the 
differential conductivity becomes a dissymmetrical function of the voltage. 

To estimate A<Jd{V), we observe that it is zero when V = 0, because n (p = pp) = 1/2 as it 
should be and it follows from Eq. (2.1) as well. It is seen from Eq. (10.4) that A<7d(V) is an even 
function of both [z — //) and V. Therefore we can assume that at low values of the voltage V the 
dissymmetrical part behaves as A<Jd(V) oc V 2 . Then, the natural scale to measure the voltage is 2T, 
as it is seen from Eq. (10.2). In fact, the dissymmetrical part is to be proportional to (pj — Pi)/pF- 
As a result, we obtain 

Here, So(r) is the temperature independent part of the entropy (see Eq. (9.1)), c is a constant, which 
is expected to be of the order of a unit. This constant can be evaluated by using analytical solvable 
models. For example, calculations of c within a simple model, when the Landau functional E[n(p)} 
is of the form 

give c ~ 1/2. It follows from Eq. (10.5), that when V ~ 2T and FC occupies a noticeable part of 
the Fermi volume, (pf —Pi)/pF — 1 5 the dissymmetrical part becomes comparable with differential 
tunnelling conductivity, Aad{V) ~ Vd(V). 

The dissymmetrical behavior of the tunnelling conductivity can be observed in measurements 
on both high-T c metals in their normal state and the heavy fermion metals, for example, such as 
YbRh 2 (Sio.95Geo.o5)2 or YbRh 2 Si2 which are expected to have undergone FCQPT. In the case of HF 
metals, upon the application of magnetic field B the effective mass is to diverge as given by Eq. (5.5). 
Here B c0 is the critical magnetic field which drives the HF metal to its magnetic field tuned quantum 
critical point. The value of the critica 
observations collected on these metals 
applying magnetic field B c0 at temperatures T*(B) < T < Tf. In the case of these metals, Tf is of 
the order of few Kelvin. We note that at sufficiently low temperatures, the application of magnetic 
field B > B c q leads to the restoration of the LFL behavior with M*(B) gi ven by Eq. (5.5). As a 



ex pone nt a = —1/2 is in good agreement with experimental 
lolllO^ . The measurements of A<Jd{V) have to be carried out 



126|. 



result, the dissymmetrical behavior of the tunnelling conductivity vanishes 

The dissymmetrical differential conductivity A<Jd{V) can also be observed when both the high-T c 
metal and the HF metal in question go from normal to superconducting. The reason is that n (p) 
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is again responsible for the dissymmetrical part of <7d{V). As we have seen in Section 3, this no(p) 
is not appreciably disturbed by the pairing interaction which is relatively weak as compared to the 
Landau interaction forming the distribution function n (p). In the case of superconductivity, we have 
to take into account that the ratio, 

^1= m nor) 

N(0) V£ 2 - A 2 ' 1 ' } 

comes into the play because the density of states N S (E) of the superconducting metal is zero in the 
gap, that is when \E\ < |A|. Here E is the quasiparticle energy, while the normal state quasiparticle 



energy is e — fi = y/E 2 — A 2 . Now we can adjust Eq. (10.4) for the case of superconducting HF 
metal, multiplying the right hand side of Eq. (10.4) by N s /N(0) and replacing the quasiparticle 
energy z — /i by \/E 2 — A 2 with E being represented by the voltage V . As a result, Eq. (10.5) can 
be presented in the following form 

, 2 r- 



Aa d (V) 



V 2 - A 2 



Pf-Pi 



1 



"A" 


2 


v- 


y. 




A. 



5o(r) (10.* 



Xpc 



|A|VV 2 - A 2 p F 

It is seen from Eqs. (10.5) and (10.8) that, as in the case described by Eq. (9.7) when the abrupt 
change in the Hall coefficient is defined by So(r), the dissymmetrical part of the differential tunnelling 
conductivity is also proportional to the Sq(t) vanishing at r —>■ 0. Note that the scale 2T entering 
Eq. (10.5) is replaced by the scale A in Eq. (10.8). In the same way, as Eq. (10.5) is valid up to 
V ~ 2T, Eq. (10.8) is valid up to V ~ 2|A|. It is seen from Eq. (10.8) that the dissymmetrical part 
of the differential tunnelling conductivity becomes as large as the differential tunnelling conductivity 
at V ~ 2|A| provided that FC occupies a large part of the Fermi volume, (pf — Pi)/pF — 1- I n the 
case of a (i-wave gap, the right hand side of Eq. (10.8) has to be integrated over the gap distribution. 
As a result, AaJ V) is expected to be finite even at V = A 1; where A± is the maximum value of the 



c?-wave gap 
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The presence of an electronic inhomogeneity in Bi 2 Sr 2 CaCu20a_|_ a _was recently discovered in ob- 
servations using scanning tunnelling microscopy and spectroscopy |l28j ] . This inhomogeneity is man- 
ifested as spatial variations in the local density of states (LDOS) spectrum, in the low-energy spectral 
weight, and in the magnitude of the superconducting energy gap. The inhomogeneity observed in 
the integrated LDOS is not induced by impurities, but rather is intrinsic in nature. The observations 
allow ed to relate the magnitude of the integrated LDOS to the local oxygen doping concentration 



128j. Spatial variation of the tunnelling differential conductance spectrum are shown in Fig. 19. The 



curves 1 and 2 can be considered as corresponding to the normal state of high-T c superconductor. 

64 



-100 -50 50 100 

Bias (mV) 



FIG. 19: Spatial variation of the tunnelling differential conductance spectra measured on Eh^S^CaC^Og+x- 
Curves 1 and 2 are taken at positions where the integrated LDOS is very small. The low differential 
conductance and the absence of a superconducting gap are indicative of insulating behavior. Curve 3 is for 
a large gap of 65 meV, with low coherence peaks. The integrated value of the LDOS at the position for 
curve 3 is small but larger than those in curves 1 and 2. Curve 4 is for a gap of 40 meV, which is close to 
the mean value of the gap distribution. Curve 5, taken at the position with the highest integrated LDOS, 
is for the smallest gap of 25 meV with two very sharp coherence peaks 

The other curves can be viewed as corresponding to high-T c superconductors with different oxygen 
doping concentrations. It is seen from Fig. 19 that the tunnelling differential conductivity is strongly 
dissymmetrical in both the normal state and the superconducting state of the E^S^CaC^Og+x com- 
pound. Therefore, we can conclude that the dissymmetrical tunnelling described by Eqs. (10.5) and 
(10.8) is in good qualitative agreement with the experimental facts displayed in Fig. 19. 

The dissymmetrical tunnelling conductivity given by Eqs. (10.5) and (10.8) can be observed in 
measurements on the heavy fermion metal CeCoIns in its superconducting state and its NFL normal 
state. As it was discussed in Subsection 9.2, this metal is expected to have undergone FCQPT. 

11. SUMMARY AND CONCLUSION 

Through out this paper we have discussed the manifestations of the fermion condensation, which 
can be compared to the Bose-Einstein condensation. A number of experimental evidences have been 
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presented that are supportive to the idea of the existence of FC in different natural and artificial 
substances. We have demonstrated also that numerous experimental facts collected in a whole variety 
of materials, belonging to the high-T c superconductors, heavy fermion metals and strongly correlated 
2D structures, can be explained within the framework of the theory based on FCQPT. 

We have shown that FCQPT separates the regions of LFL and strongly correlated liquids. Be- 
yond the FCQPT point the quasiparticle system is divided into two subsystems, one containing 
normal quasiparticles, while the other being occupied by FC localized at the Fermi level. In the 
superconducting state the quasiparticle dispersion e(p) in systems with FC can be approximated by 
two straight lines, characterized by effective masses M* FC and M£, and intersecting near the binding 
energy E which is of the order of the superconducting gap. The same quasiparticle picture and the 
energy scale E persist in the normal state. 

We have demonstrated that fermion systems with FC have features of "quantum protectorate" 
and shown that the theory of high-T c superconductivity, based on FCQPT and on the conventional 
theory of superconductivity, permits the description of high values of T c and of the maximum value 
of the gap Ai, which may be as big as Ai ~ O.Ief or even larger. We have also traced the tran- 
sition from conventional superconductors to high-T c ones. We have shown by a simple, although 
self-consistent analysis that both the pseudogap state and the general features of the shape of the 
critical temperature T c (x) as a function of the number density x of the mobile carriers in the high-T c 
compounds can be understood within the framework of the theory. 

We have also shown that striking experimental results on the magnetic-field induced LFL in high- 
Tc metals, which unveil the nature of the high-T c superconductivity, suggest that FCQPT and the 
emergence of the novel quasiparticles with effective mass strongly depending on the magnetic field 
and temperature and resembling the Landau quasiparticles are qualities intrinsic to the electronic 
system of the high-T c superconductors. 

We have provided explanations of the experimental data on the divergence of the effective mass 
in 2D electron liquid and in 2D 3 He, as well as shown that above the critical point of FCQPT the 
system exhibits the LFL behavior. The behavior of the heavy electron liquid approaching FCQPT 
form the disordered phase can be viewed as the highly correlated one because the effective mass is 
very large and strongly depends on the density, temperature and magnetic fields. 

At different temperatures, the behavior in magnetic fields of the highly correlated electron liq- 
uid approaching FCQPT from the disordered phase has been considered. We have shown that at 
sufficiently high temperatures T\(B) < T, the effective mass starts to depend on T, M* oc T -1 / 2 . 



66 



This T" 1 / 2 dependence of the effective mass at elevated temperatures leads to the non-Fermi liquid 
behavior of the resistivity, p(T) oc T. The application of magnetic field B restores the common 
T 2 behavior of the resistivity. If the magnetic field (B — B c0 ) decreases to zero, the effective mass 
M* diverges as M* oc (B — B c0 )~ 2 / 3 . At finite magnetic fields, the regime NFL is restored at some 
temperature Ti(B) oc (B — B c0 ) with M*{T) oc T~ 2 / 3 . We have demonstrated that this B — T phase 
diagram has a strong impact on MR of the highly correlated electron liquid. At fixed B, MR as a 
function of the temperature exhibits a transition from the negative values of MR at T — > to the 
positive values at T oc (B — B c0 ). While at low temperatures and elevated magnetic fields, MR goes 
from positive to negative. This behavior was observed in the heavy fermion metals. 

We have demonstrated that the strongly correlated electron liquid with FC, which exhibits strong 
deviations from the LFL behavior down to lowest temperatures, can be driven into LFL by applying 
a magnetic field B. If the magnetic field (B — B cQ ) decreases to zero, the effective mass M* diverges 
as M* oc 1/y/B — B c o and the Neel temperature T/v of the AF phase transition tends to zero, 
T N (B — > B c0 ) — > 0. The NFL regime is restored at some temperature T*(B) oc \[B — B c0 . In that 
case and at T — > and B = B c0 , the Griineisen ratio as a function of temperature T diverges. While 
the entropy S(T) possesses the specific low temperature behavior, S(T) oc S + aVT + bT with So, a 
and b are temperature independent constants. We have shown that the obtained T — B phase diagram 
is in good agreement with the experimental T—B phase diagram obtained in measurements on the HF 
metals YbRh 2 Si2 and YbRh 2 (Sio.95Geo.o5)2- We have also demonstrated that the abrupt jump in the 
Hall coefficient Rh(B — > B c0 , T — > 0) is determined by the presence of FC. We have observed that the 
second order AF phase transition changes to the first order one below T crit <T N (B = 0) making the 
corresponding quantum and thermal critical fluctuations vanish at the jump. Therefore, the abrupt 
jump and the divergence of the effective mass taking place at T/v — > are defined by the behavior of 
quasiparticles rather than by the corresponding thermal and quantum critical fluctuations. 

We have predicted that the differential tunnelling conductivity between a metallic point and an 
ordinary metal, which is commonly symmetric as a function of the voltage, becomes noticeably 
dissymmetrical when the ordinary metal is replaced by a HF metal, the electronic system of which 
has undergone FCQPT. This dissymmetry can be observed when the HF metal is both normal 
and superconducting. We have also discussed possible experiments to study the dissymmetry in 
measurements on the HF metals. In the case of the high-T c superconductors, our consideration is in 
good agreement with available data. 

In conclusion, we have shown that in contrast to covnetional quantum phase transitions, whose 



67 



physics is dominated by thermal and quantum fluctuations and characterized by the absence of 
quasiparticles, the physics of Fermi systems and the heavy electron liquid near FCQPT or undergone 
FCQPT is determined by quasiparticles resembling the Landau quasiparticles. We have shown that 
the Landau paradigm based on the notion of quasiparticles and order parameters is still applicable 
when considering the low temperature properties of the heavy electron liquid, whose understand- 
ing has been problematic largely because of the absence of theoretical guidance. In contrast with 
the conventional Landau quasiparticles, the effective mass of the considered quasiparticles strongly 
depends on the temperature, applied magnetic fields, the number density x, pressure, etc. These 
quasiparticles and the order parameter are well defined and capable of describing both the LFL and 
the NFL behaviors of both the high-T c superconductors and the HF metals and their universal ther- 
modynamic properties down to the lowest temperatures. This system of quasiparticles determines 
the recovery of the LFL behavior under applied magnetic fields and preserves the Kadowaki- Woods 
ratio. Thus, we obtain a unique possibility to control the essence of HF systems and accordingly of 
the HF metals by magnetic fields in a wide range of temperatures. 

Finally, our general consideration suggests that FCQPT and the emergence of novel quasiparticles 
at QCP and behind QCP and resembling the Landau quasiparticles are qualities intrinsic to strongly 
correlated substances, while FCQPT can be viewed as the universal cause of the non-Fermi liquid 
behavior observed in different metals and liquids. 
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